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ABSTRACT 

We present an analysis of the molecular and atomic gas emission in the rest-frame 
far-infrared and sub-millimetre, from the lensed z = 2.3 sub-millimetre galaxy 
SMM J2135— 0102. We obtain very high signal-to-noise detections of 11 transitions 
from three species and limits on a further 20 transitions from nine species. We use the 
12 CO, [Ci] and HCN line strengths to investigate the gas mass, kinematic structure and 
interstellar medium (ISM) chemistry, and find strong evidence for a two-phase medium 
within this high-redshift starburst galaxy, comprising a hot, dense, luminous compo- 
nent and an underlying extended cool, low-excitation massive component. Employing 
a suite of photo-dissociation region models we show that on average the molecular gas 
is exposed to a UV radiation field that is ~ 1000 x more intense than the Milky Way, 
with star-forming regions having a characteristic density of n ~ 10 4 cm~ 3 . Thus, the 
average ISM density and far-UV radiation field intensity are similar to those found in 
local ULIRGs. These densities and radiation fields are similar to those found in the 
central regions of typical starburst galaxies, even though the star formation rate is far 
higher in this system. The 12 CO spectral line energy distribution and line profiles give 
strong evidence that the system comprises multiple kinematic components with dif- 
ferent conditions, including temperature, and line ratios suggestive of high cosmic ray 
flux within clouds, likely as a result of high star formation density. We find tentative 
evidence of a factor ~ 4 temperature range within the system. We expect that such 
internal structures are common in high-redshift ULIRGs but are missed due to the 
poor signal-to-noise of typical observations. We show that, when integrated over the 
galaxy, the gas and star-formation surface densities appear to follow the Kennicutt- 
Schmidt relation, although by comparing our data to high-resolution sub- mm imaging, 
our data suggest that this relation breaks down on scales of <100pc. By virtue of the 
lens amplification, these observations uncover a wealth of information on the star for- 
mation and ISM at z ~ 2.3 at a level of detail that has only recently become possible 
at z < 0.1, and show the potential physical properties that will be studied in unlensed 
galaxies when the Atacama Large Millimeter Array is in full operation. 

Key words: galaxies: active — galaxies: evolution — galaxies: high-redshift — galax- 
ies: starburst — sub-millimetre 
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1 INTRODUCTION 

The properties of the interstellar medium (ISM) play a criti- 
cal role in the evolution of galaxies as it includes the raw ma- 
terial from which stars form. The main sources of heating of 
the atomic and molecular ISM are the stellar radiation field, 
cosmic rays and turbulence, and as such the ISM exhibits a 
considerable range of properties, spanning a wide range in 



density 



10" 



cm" d and temperature ~ 10-1000 K. The 



with star formation rates > 10" 



19981. 1 19991; iGreve et al. 


12005 


Hailev-Dunsheath et al. 


2010a 



thermodynamic state of the gas is dictated by the balance of 
this heating with cooling from atomic and molecular species 
(e.g. [Cn], [Cl], 12 CO and [Ol]) and this balance may dictate 
the Jeans mass and initial mass fu nction (IMF) and there - 
fore the efficiency of star formation l|Hocuk fc SpaansfeoioT ) . 
Understanding the balance of heating and cooling within 
the ISM is thus fundamental to understanding the detailed 
physics of the star formation which drives the formation and 
evolution of galaxies. 

An essential element of our understanding of distant 
galaxies comes therefore from observations of emission from 
interstellar molecular and atomic gas. To date, the ma- 
jority of the detections of 12 CO and atomic fine struc- 
ture emission from high-redshift galaxies have been in pow- 
erful QSOs and intrinsically luminous galax ies, typically 

M fl yr" 1 (e.g. iFraver et all 
Tacconi et all 120061 . |200S| ; 
Walter et al.Tl2009l ). How- 
ever, the detailed physical properties of the vigorous star- 
bursts within high-redshift starburst galaxies are still un- 
known. Analogy to local ultra- luminous infrared galaxies 
(ULIRGs) would argue for compact starbursts triggered 
by mergers, although secular bursts in massive gas disks 
are also feasible (e.g. iDave et all |2010| ). These different 
scenarios imply different properties for the ISM and so 
may be tested using observations of the molecular and 
atomic emission. Critically, the first full 12 CO spectral line 
energy distributions (SLEDs), including high-J uppcr tran- 
sitions, have recently been completed for local ULIRGs 
l|Papadopoulos et all2Qlol ; lvan der Werf et al J2010h provid- 
ing an empirical benchmark for comparison to high-redshift 
sources, to search for similarities in the thermodynamics of 
their ISM. 

Much of the atomic and molecular gas emission from 
galaxies arises from photo-dissociation regions (PDRs), the 
surfaces of molecular clouds where the heating and chem- 
istry is dominated by far-ultraviolet (far-UV) photons from 
stellar sources (Hollenbach & Tielens 1997). PDRs dominate 
the infrared and sub-millimetre emission line spectra of star 
formation regions and galaxies as a whole. The majority of 
their cooling occurs via atomic fine structure lines o f [Ol], 
[Cn] and the rotational 12 CO lines l|Kramer et alll2004 ). 

Prior to full science operations of ALMA, the most 
promising route to gaining high signal-to-noise, sub- 
millimetre spectroscopy of star-forming galaxies at high red- 



tational lensing (e.g. 


Smail et all 20021; Weifi et al. 120051 ?L 


20071; iMaiolino et all 


20091; lHailev-Dunsheath et all l2010al; 


Ivison et all 120101). In particular, measuring the chemistry 



of the ISM through emission line flux ratios in strongly- 
lensed galaxies is particularly advantageous since gravita- 
tional lensing is achromatic (in the absence of differential 
amplification) and so the apparent line ratios (and the im- 



plied chemistry) are unaffected. As such, intrinsically faint 
emission lines which are sensitive probes of the ISM chem- 
istry (such as 13 CO or HCN) can be measured, if suitably 
bright targets can be identified . 

Recently. ISwinbank et all (|2010T ) reported the discovery 
of a strongly-lensed sub- millimetre galaxy SMM J2135— 0102 
(hereafter SMMJ2135; R.A./Dec: 213511.6, -0102 52; 
J2000) behind a massive galaxy cluster at z = 0.325. The 
redshift for the galaxy, z = 2.3259, was derived through 
blind detection of 12 CO(1-0) using Zpectrometer on Green 
Bank Telescope (GBT). Based on a detailed lensing model, 
the 870-/im emission appears to be amplified by a factor 
32.5 ± 4.5x, implying an unle nsed flux density o f ~ 3mJy 
(Lboi = 2.3 ± 0.4 x 10 12 L Q ; llvison et ai1l2010f ). which is 
close to the confusion limit of existing sub-millimetre sur- 
veys (|Blain et al.ll2002l ). 

This uniquely bright source therefore provides an op- 
portunity to investigate the detailed properties of the ISM 
within a galaxy which is representative of the high-redshift 
starburst population. In this paper we analyse follow-up ob- 
servations of the 12 CO ladder (up to J up per = 10), as well 
as other dense gas tracers and more complex molecules in 
the rest-frame frequency range 80-1100 GHz. To put the 
current molecular/atomic line data obtained for the global 
line emission of this distant starburst in perspective, we 
note that they are rivaled only by those obtai ned for local 
ULIRGs, e.g. Arp220, Mrk231 and NGC6240 (|Greve et all 
120091 ; Ivan der Werf et all lioich . which are some ~ 250 x 
closer to us. 

In Sj2]we describe our observations of the molecular and 
atomic emission from SMM J2135. Our analysis is described 
in f|3] where we first examine the integrated properties of 
the system, deriving and comparing gas masses calculated 
from 12 CO, [Ci] and HCN. We produce a 12 CO SLED for 
this galaxy and use this to investigate excitation conditions 
within its ISM. We use PDR models to investigate the char- 
acteristic density of the ISM and the incident far-UV in- 
tensity. We then decompose the spectra into three distinct 
kinematic components, and measure the 12 CO SLED of each 
in order to search for excitation structure within the system. 
We give our conclusions in ^4] Throughout the paper we use 
a ACDM cosmology with Hp = 72kms~ : *- Mpc" 1 , ft m = 0.27 
and Qa = 1 — (|Spergel et all l2003h . Unless otherwise 
stated, a lensing amplification correction of a factor of 32.5 
has already been applied to all quoted luminosities. 



2 OBSERVATIONS 

An essential requirement for the study of the molecular and 
atomic emission lines from SMM J2135 is a precise redshift 
for the gas reservoir in this system. This was obtained soon 
after the discovery of t he source using Zpectrometer on GBT 
jSwinbank et alj201(j see j|2Tj| . With this accurate systemic 
redshift, it was then possible to precisely tune to the ex- 
pected frequencies of molecular and atomic emission lines. 



2.1 GBT Zpectrometer Observations 

Details of the observa tions with Zpectrometer are given in 
ISwinbank et all (|2010h and Table 1. Briefly, Zpectrometer is 
a wide-band spectrometer optimised for 12 CO(1-0) emission 



© 0000 RAS, MNRAS 000, 000-000 



The ISM in a z ~ 2.3 Star-Forming Galaxy 3 



Table 1. 

Log of IRAM 30 m and GBT Observations 



Band iv obs Emission lines ti nt 

(GHz) (ks) 



Ka 


25.6- 


-36.1 


12 CO(1-0), CS(2-1), CN(l-O) 










HNC(l-O), 13 CO(1-0) 










HCO+(1-0) 


18.0 


E0 


79.00- 


-83.00 


HCN(3-2), HCO+(3-2), 










HNC(3-2) 


28.8 




96.75- 


100.75 


13 CO(3-2), H 2 325.141 


18.0 




101.97- 


-105.97 


12 CO(3-2), CN(3-2), CS(7-6) 


18.0 


El 


131.00- 


-135.00 


13 CO(4-3), HCO+(5-4), 










HCN(5-4), CS(9-8) 


14.4 




135.92- 


-139.92 


12 CO(4-3), CN(4-3), HNC(5-4) 


10.8 




146.00- 


-150.00 


[Ci](l-0), CS(10-9), 2 487.3 


7.2 




171.26-174 


12 CO(5-4) 


10.8 


E2 


205.90- 


-209.90 


12 CO(6-5) 


8.4 




224.11- 


-228.11 


H 2 731.681 


3.6 




240.53- 


-244.53 


12 CO(7-6), [Ci](2-1) 


12.6 


E3 


275.16- 


-279.16 


12 CO(8-7) 


9.6 




309.77- 


-313.77 


12 CO(9-8) 


7.2 



The log of the IRAM 30m and GBT Zpectrometer observations 
giving the frequency ranges of Zpectrometer Ka band and the 
IRAM/EMIR receiver backends (E0-E3), the emission lines ob- 
served and integration time for each set-up. 

Note: These set-ups have beam sizes of Ka:- 15"-23", E0:~ 30", 
El:~ 17", E2:~ ll"and E3:~ 8.5"respectively. 

line searches between z = 2.2-3.5 (with ~ 150 km s~ reso- 
lution ) us ing the GBT's Ka- band receiver (see lHarris et al.l 
2007, and lHarris et all [2010 for more detail on the instru- 
ment and observing mode). Observations of SMMJ2135 
were conducted in two equal shifts on 2009 May 19 and 
May 27 for a total integration time of 5 hours. Data reduc- 
tion was carried out with the standard Zpectrometer GBT 
reduction scripts. The final spectrum reaches an rms noise 
level of a = 0.50 mjy beam -1 . From the spectrum, shown in 
Fig.[TJ we determine the heliocentric redshift as z — 2.32591 
from the 12 CO(1-0) emission at 34.648 GHz. The observed 
velocity-integrated flux in 12 CO(1-0) is given in Table 2. 

2.2 IRAM PdBI & 30m Observations 

We used the six-element IRAM Plateau de Bure Interfer- 
ometer (PdBI) to observe the redshifted 12 CO (3-2) and 
12 CO(4-3) emission lines and the continuum at 104 GHz and 
139 GHz respectively. Observations were made in D configu- 
ration in Directors Discretionary Time (DDT) on 2009 May 
29 and 31 with good atmospheric phase stability and pre- 
cipitable water vapour (seeing = 0.6-1.6", pwv = 5-15mm). 
The receivers were tuned to the systemic redshift deter- 
mined from the GBT 12 CO (1-0) spectrum. We observed 
SMM J2135 for total on-source exposure times of 4hrs and 
2hrs for 12 CO(3-2) and 12 CO(4-3) respectively. The cor- 
relator was adjusted to a frequency resolution of 2.5 MHz, 
yielding 980-MHz coverage. The overall flux scale for each 
observing epoch was set from observations of MWC 349, 
with additional observations of 2134+004 for phase and am- 
plitude calibrations. The data were calibrated, mapped and 



Table 2. 
Emission Line Properties 



Species 




Flux ' 6 


7 / 

L 




(GHz) 


{ T^r tin Q \ 


fin 8 K km s~ 1 nr 2 l 


12 CO(1-0) 


115.2712 


2.16 ±0.11 


173 ±9 


12 CO(3-2) 


345.7959 


13.20 ±0.10 


117.6 ±0.9 


12 CO(4-3) 


461.0408 


17.3 ± 1.2 


87 ±6 


12 CO(5-4) 


576.2679 


18.7 ± 0.8 


60 ± 3 




oyi.4 / oi 


Zl.O ± 1.1 


join 
4o ± Z 


12 CO(7-6) 


806.6518 


12.6 ± 0.6 


21 ± 1 


12 CO(8-7) 


921.7997 


8.8 ± 0.5 


11 ± 1 


12 CO(9-8) 


1036.9124 


3.6 ±1.7 


4±2 


12 CO(10-9) 


1151.9855 


< 1.1 


< 0.9 


13 CO(1-0) 


110.2013 


< 0.3 


< 29 


13 CO(3-2) 


330.5880 


< 1.8 


< 18 


13 CO(4-3) 


440.7651 


< 0.3 


< 1.4 


[Ci](l-0) 


492.1607 


16.0 ±0.5 


71 ±2 


[Cl](2-1) 


809.3435 


16.2 ± 0.6 


26 ± 1 


HCN(l-O) 


88.6300 


< 0.3 


< 45 


HCN(3-2) 


265.8900 


1.2 ±0.4 


18 ±6 


HCN(5-4) 


443.1200 


< 0.3 


< 1.4 


HNC(l-O) 


90.6600 


< 0.3 


< 43 


HNC(3-2) 


271.9800 


< 1.2 


< 17 


HGO+(1-0) 


89.190 


< 0.3 


< 44 


HCO+(3-2) 


267.560 


< 1.2 


< 17 


HCO+(5-4) 


445.90 


< 0.3 


< 1.3 


CN(l-O) 


113.320 


< 0.3 


< 27 


CN(3-2) 


339.450 


< 1.6 


< 14 


CS(2-1) 


97.981 


< 0.3 


< 37 


CS(7-6) 


342.883 


< 1.6 


< 14 


CS(9-8) 


428.875 


< 0.3 


< 14 


CS(10-9) 


489.751 


< 1.6 


< 7 


o 2 


487.246 


< 1.6 


< 7 


H 2 


325.141 


< 1.8 


< 18 


H 2 


731.681 


< 4.8 


< 9.5 


[Cn] c 


1910 


850 ± 180 


2500 ± 500 


[OI] d 


2065.40 


620 ± 200 


155 ± 50 



a We quote 3-c limits for all lines which are not formally detected. 
'Uncertainties on fluxes denote measurement errors and do not in- 
clude the calibration uncertainties, which we estimate as ~ 5% for 
30-200 GHz, ~ 10% for 200- 300 GHz and ~ 15% for > 300 GHz. 
c The [Cul flux is ta ken from llvison et al.l J201dh . 
" Ilvison et al ] J- 1 " III report a ~ 3<r from [OI]146/im and we have 
used the flux of this feature in our analysis, although we note that 
removing this line from our analysis changes our results negligi- 
bly. We also note that the since the [Cn] and [Oi] 146/im lines are 
of comparable strength, this would imply that the [Oi] 63/xm line 
would be much brighter than the [Cn] line thus making this an 
interesting source for Herschel/PACS follow-up. 

analyzed using the GILDAS software package. Inspection of 
the velocity datacubes shows very good detections of both 
12 CO(3-2) and 12 CO(4-3) emission lines (S/N ~ 300 in each) 
at the position of SMM J2135 (see Fig. [1} and we give line 
fluxes in Table 2. 

To constrain the high-J upper 12 CO emission (and search 
for other emission lines such as [Ci]) we used the Eight 
Mixer Receiver (EMIR) multi-band hetrodyne receiver at 
the IRAM 30-m telescope. Observations were made on 2009 
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June 29-30, and 2010 February 5-8 and April 3-6, in good 
to excellent conditions, typically with ~ 2-6 mm pwv and 
;£ 1 mm pwv for some of the April observations. Data were 
recorded using the E0-E3 receivers with 4 GHz of instan- 
taneous, dual-polarisation band-width covering frequency 
ranges from ~ 70-310 GHz (see Table 1 for the details of the 
setups). For each observation, we used eight 1-GHz band- 
width units of the Wide-band Line Multiple Autocorrelator 
(WILMA) to cover 4 GHz in both polarizations. WILMA 
provides a spectral resolution of 2 MHz which corresponds to 
5-7 km s -1 in the 3-mm band. The observations were carried 
out in wobbler-switching mode, with a switching frequency 
of 1 Hz and an azimuthal wobbler throw of 90" . Pointing was 
checked frequently on either the nearby quasar J2134+004 
or Venus and was found to be stable to within 3". Cali- 
bration was carried out every 12 mins using the standard 
hot/cold- load absorber, and the flux calibration was carried 
out using the point source conversion between temperature 
and flux as measured from celestial objects. We note that, 
even at the highest frequency (i.e. 345 GHz) the beam is 8.5" 
and so the lensed gal axy is likely to be smal ler than the beam 
at all wavelengths (|Swinbank et al.l l20ld). The data were 
first processed with the CLASS software l|Pety||2005h . and 
then using custom idl routines. We omitted scans with dis- 
torted baselines and subtracted only linear baselines from in- 
dividual spectra. For the 2010 April observations of 12 CO(5- 
4), an upward correction was made to the flux of 10% since 
the line extends into the wings of an atmospheric water band 
(although for these observations the pwv g 1mm). In total, 
we observed twelve different set-ups, each for ~ 2-5 hrs (Ta- 
ble 1); the spectra of the detections are shown in Fig. [T] 
These observations provide detections or limits on 31 tran- 
sitions listed in Table 2. 



2.3 APEX SHFI Observations 

Observations of the redshifted 12 CO (7-6) emission line were 
also carried out in DDT using the Swedish Hetrodyne Fa- 
cility Instrument (SHFI) between 2009 July 10 and July 
20 (ESO programme ID 283.A-5014). SHFI consists of four 
wide-band heterody ne receiver channels for 230-1300 GHz 
l|Vassilev et al.l |200Sj). We used APEX-1 tuned at 242.4 GHz 
to search for the 12 CO(7-6) emission from the galaxy, and 
obtained a total integration time of ~ 5.5 hrs in excellent 
conditions (0.4-0.6 mm pwv). The data reduction was car- 
ried out using the CLASS software, omitting scans with poor 
baselines. To create the final spectrum we binned the data 
onto a velocity scale of 50 kms -1 and averaged it together 
with 4 hours of observations from IRAM/EMIR, which has 
similar signal-to-noise; we show this spectrum in Fig. [TJ 



in Fig. [TJ The resulting combined spectrum reaches an rms 
of 2mJy per 100 km s -1 channel and although 12 CO (10- 
9) is not detected, we place a 3-cr limit on the flux of 
< 0.2 Jykms" 1 . 



3 ANALYSIS AND DISCUSSION 

Our observations have detected 11 individual transitions 
arising from three molecular or atomic species in SMM J2135 
and place upper limits on the line fluxes of a further 20 
transitions arising from seven other species. We list these in 
Table 2 with the fluxes quoted with their respective mea- 
surement uncertainties and we show the spectra for all de- 
tections in Fig. [1] We estimate additional calibration uncer- 
tainties of 5%, 10% and 15% for those lines at 30-200 GHz, 
200-300 GHz and > 300 GHz respectively. 

With the very high signal-to-noise of our data (partic- 
ularly the 12 CO lines), we clearly identify multiple velocity 
components within the spectra (at approximately —180, 
and +200 kms -1 relative to the 12 CO(1-0) flux-weighted 
mean velocity), implying multiple physical components in 
the source. We believe that the detection of this structure 
in the 12 CO spectra is not due to this galaxy's being un- 
usual, but instead reflects the high signal-to-noise of our 
observations, which is uncommon for observations of other 
(even local) galaxies. We discuss these multiple components 
in £13 . 3 1 modelling the system as three components which we 
label X, Y and Z. However, so that fair comparisons can be 
made to other galaxies (at low and high redshift), we begin 
(in i]3.1.1l 13.1.21 and 13. 2p by discussing the integrated prop- 
erties of the system to see what can be learnt about the bulk 
properties of its ISM. 

As Fig. [1] shows, the 12 CO and [Cl] emission lines are 
broad, with a typical FWHM of ~ 500 km s -1 . Since there 
is significant kinematic structure within the 12 CO and [Cl] 
emission, we first define the width of the emission lines 
by constructing a composite 12 CO spectrum by combining 
the spectra from 12 CO (1-0) to 12 CO (8-7) (normalised by 
peak flux). This composite has a full width at zero inten- 
sity of 900 kms -1 , with a range from —350 to 550 kms" 
(Fig. [TJ), although we note that there may be faint emis- 
sion extending to ±1000km s . We use the velocity range, 
—350 to 550 kms -1 , to measure the fluxes of the 12 CO lines, 
[Cl]( 3 Pi ->■ 3 P ) and [Cl]( 3 P 2 3 Pi), while for HCN(3-2) 
we use a range of —100 to 550 kms -1 . To determine the 
errors on the fluxes we measure the variance in the spec- 
tra away from the emission line, although in most cases the 
dominant error is from the calibration; these errors are used 
in Fig. [2] (see also Table 2). 



2.4 SMA Observations 

We mapped the 870-^tm continuum emission from 
SMM J2135 using the Submillimeter Array ( SMA) in a num- 
ber o f different array configurations (see Swinban k et al.l 
|2010T ). During four of these tracks (in the sub-compact, com- 
pact, extended and very-extended configurations) we tuned 
the upper side band of the receiver to 346.36699 GHz to 
search f or 12 CO (10-9). Detail s of the data reduction are 
given m ISwinbank et all (|2010h and we show the spectrum 



3.1 Integrated Properties 

3.1.1 Gas Mass and Dust SED 

We begin by estimating the total molecular gas mass 
of the galaxy from 12 CO. Low-J uppcr transitions of 
12 CO are commonly used to obtain estimates of the 
total gas mass of a galaxy using a conversion factor, 
a = Mh 2 / L'i2 CO (i-o) (K kms" 1 pc 2 ) -1 which converts 
the 12 CO line luminosity to total gas mass, (where Mh 2 is 
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Figure 1. Molecular and atomic emission from the lensed sub-millimetre galaxy SMM J2135. The top two rows show spectra of the 1 CO 



emission arising from J u 



1 up to J u p 



10. These are followed by the spectra of the [Cl] fine-structure lines and our detection 



of HCN(3— 2). In all cases, the emission profiles show multiple velocity components whose intensity appears to vary between transitions. 
To better constrain the kinematic structure of the lines we derive an average 12 CO spectrum (which does not include 12 CO(9-8) or 
12 CO(10-9)) and we overplot on this the resulting three-component parametric model as described in £13,31 We fit this three-component 
kinematic model to the various lines, allowing the intensities of the components to vary between lines, and overlay this on all the 12 CO 
emission lines making up the average, clearly highlighting the evidence of different excitation in the different velocity components. The 
spectra have been binned into channels of width 30-100 km s _1 . The spectra for 12 CO(3-2), 12 CO(4-3) and 12 CO(6-5) are the combined 
data from our PdBI and EMIR observations, the 12 CO(7-6) spectrum is a combination of EMIR and SHFI data. The 12 CO(1-0) line 
was measured by GBT and the 12 CO(10-9) line by SMA, all other observations were taken with IRAM 30m and PdBI. The axes are 
flux (mjy) and frequency (GHz) (with the velocity shown on the upper horizontal axes) in all cases except for the average spectrum 
which has been plotted as kms" 1 and flux normalised to the peak value. Note that 12 CO(7-6) and [Cl](2— 1) lines abut each other in 
the spectrum. 



def ined to include the mass of Heliu m such that Mh, = M gaa ; 
see lSolomon fc Vanden BoutJl2005l for a review). 

There is considerable uncertainty about the correct 
value of q to adopt in high-redshift galaxies (e. g. iBaker et al.l 
12004 ICoppin et all 120071 : iTacconi et afl|20Qg| and Ivison et 
al. 2010c) and similar uncertainty about the conversion fac- 
tor to transform high-Ju PP cr CO transitions to the equiva- 
lent CO (1-0) fluxes (if required). However, we can estimate 
a minimum mass (and so a lower limit on a) assuming the 
1 CO (1-0) emission is optically thin (which i s unlikely to 
be the case) and has solar abundance, following llvison et al.l 
|20ld ): 



M(H 2 



n.„.s 1 9i e -^/n ( m 



J CO(1-0) 



Z 

12 CO/H 2 



io- 



-1 



)-J(T b9 ) \ 
J(Tk) J 



Mr: 



Kkms~ 1 pc 2 



where T a = E 1 /k B ~ 5.5 K, J(T) = T (e T °/ T - 1)"\ 

T bg = (1 + z)T C mb ~ 9.1 K (Tomb = 2.73 K at z = 0), 
2 = 2.3259, gi — 3 (the degeneracy of level n = 1), Z ~ 
2(T k /T ) and [ 12 CO/H 2 1~ 10~ 4 fo r a solar metallicity envi- 
ronment l|Brvant fc Scoviiie| [l996). For typical star- forming 
gas, Tk ~ 40-60 K, so to derive a minimum mass we use 



T k = 40K. With L' 12co(1 _ 0) = 173 ± 9 x 10 8 Kkms -1 pc 2 , 
this gives us a minimum gas mass of M gas J£ 1 x 10 Mq 
and a lower limit on a £ 0.54. 

We can then compa re this to the dynamical mas s 
adopting M dyn =5 Ra 2 /G (|Solomon fc Vanden Boutll2005l ). 
This assumes the dynamics of the CO emission trace the 
virialised potential well of the system and that the CO is 
distributed i n a sphere of constant density with r ~ 1 kpc 
(see %Q1 and lSwinbank et al.ll2010l ). Using the width of the 
12 CO(1-0) line, a ~ 200 kms -1 , we estimate a dynamical 
mass of M dyn ~ 5.3 x 10 10 M s . If this mass is dominated by 



gas, then the upper limit on a is a = M gas /L^ 



12 CO(1-0) 



< 3. 



Using this dynamical mass and the minimum gas 
mass, we place a lower limit on the gas fraction of 
M gas /Mdyn k, 0.19. This ratio is similar to the ratio for typ- 
ical starburst nuclei a nd the ratio of 30% found in M 82 
ijDevereux et al.l 1 19941 ). Compared to other high-redshift 
sub- millimetre galaxies (SMGs), u sing the 12 CO (3-2) ob- 
servations from iGreve et ail (|2005l ) and adopting a = 0.8 
and R 3 ,i = 0.58 ± 0.05 from Ivison et al. (2010c), we find 
a median gas mass r atio of 50% for luminous SMGs, whilst 
ITacconi et ail <|2008l ) derive a gas mass fraction of 0.3 for 
more typical star-forming galaxies at similar redshifts. We 
stress, however, that compared to the gas mass of 6 x 10 8 Mq 
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within the central 1.2 kpc of M82 |Young fc Scovillelll984h . 
our gas mass is nearly two orders of magnitude higher, un- 
derlining the more extreme conditions in the central regions 
of SMMJ2135. 

In order to provide a simple comparison with pre- 
vious studies, we estimate the gas mass of SMM J2135 
using a conversion factor of a ~ 0.8 which applies to 
the smoothly distributed, high-pressu re, largely molecu- 
lar ISM measured in local ULIRG s (|Downes fc Solomon! 
1 19981 ; ISolomon fc Vanden Boudl2005l . although recent local 
studies on high-J uppcr CO lines imply multi-phase ISM in 
ULIRGs as opposed to smoothly distributed) and is also the 
canonical value used for high-redshift LIRGs and ULIRGs 
jTacconi et all 120081 : IStark et alj|2008h . This value of a is 
comfortably between the upper and lower limits calculated 
above. With this assumption the 12 CO(1-0) line luminosity 
yiel ds a total gas ma ss of Mh 2 = (1.4 ±0.1) x 1O 1O M (see 
also ISwinbank et all l20ld ) . 

Taking our estimate of the total gas mass and the ex - 
pected size of ~ 1 kpc for the system l|Swinbank et aljfeoiol ). 
we derive an average column density of ~ 10 24 cm -2 , which 
is comparable to the molecular hydrogen density in Arp 220, 
averaged over a similar region |Gerin fc Phillips] 1 1998h . but 
much higher than the density in typical starburst galax- 
ies such as M 82. However, this estimate gives an average 
over the system which we know to be structured and so 
indicates that the column density is very high in parts of 
this source. The associated extinction is expected to exceed 
Ay ~ 10 3 , suggesting significant absorption even in the far- 
infrared and moreover that the emission from some of the 
more common species we see is optically thick. With this 
extinction, coupled with our high value for Go (derived in 
N3.2.ip . high infrared luminosity and potentially high cos- 
mic ray flux (see below), it is probable that momentum- 
driven outflows will result. These out flows could be either 
photon-driven (Thompson et al.ll2005h or cosmic ray-driven 
|Socrates et al.l 20081 ') and will impact both the dynamics of 



the gas in the system and the steady-state assumption in 
our PDR modelling ( H3.2.ip . 

We can now compare the star-formation efficiency in 
SMMJ2135 to that in the wider SMG and ULIRG pop- 
ulations. First we note that the far-infrared l uminosity of 
SMMJ2135 is L Fm = (2.3 ± 0.1) x 10 12 L Q (|lvison et all 



2010) which indica tes a star forma tion rate (SFR) of 
400 ± 2OM yr _1 (|Kennicuttl 1 19981 ) assuming a Salpeter 
IMF. We combine the star-formation rate and gas mass 
to estimate the st ar formation efficiency (SFE) following 
iGreve et all (|2005h . SFE= Lfir/Mh 2 resulting in SFE~ 
165 ± 7 Lq M^ 1 . This i s comfortably within the limit de- 



500 L^M 



rived by IScovilld \2004 ) of [Lfir/Mh 2 
assuming Eddington-limited accretion of gas onto OB star 
clusters. 12 CO(1-0) measurements in local ULIRGs derive 
a SFE of 180 ± 160 L0 Mq 1 , which is comparable to that in 
our source. The SFE of our source is also comparable to the 
median S FE of 210 ± 8 Lq for the sample of luminous 
SMGs of IGreve et al.l l|2005f ). after corre cting the latter us- 
ing .R3.1 = 0.58 as appropriate for SMGs (|Harris et al.ll2010l : 
Ivison et al. 2010c). 

We begin modelling this system by noting that the 870- 
/j,m SMA observations and comparison to the LABOCA flux 
reveal four bright, compact clumps embedded in a much 
more extended system with the clumps emitting ~ 80% of 



the total luminosity, and the rest of the emission emergin g 
from a more extended component l|Swinbank et al.l 12010) . 
Motivated by this structure we fit a two-component model 
to the dust SED, fixing the characteristic temperature of the 
cool (presumably more extended) component at 30 K and 
allowing the size and the characteristic temperature of the 
clumps to vary. The best fit yields a characteristic size for the 
clumps of r ~ 200 pc at a temperature of Td.warm = 57±3 K, 
while the extended component has a size of r ~ 1000 pc at 
Td.cooi = 30 K. The dust mass of both the extended compo- 
nent and the clumps is then (1.5 ± 0.2) x 10 s M©, which is 
si milar to the mass d etermined for the cool dust component 
in I Ivison et alj l|201Cf ). The inferred clump size is somewhat 
larger than the size of the dust emission regions seen in the 
highest resolution 870- /im SMA map s, where the clumps ap - 
pear to be 100-200 pc in diameter (|Swinbank et al.l l2010). 
However, we note that forcing the clump sizes down to that 
size requires an increase in the characteristic temperature 
up to ~ 80 K, and to keep the emission consistent with 
the integrated SED the clumps must then be optically thick 
at ~ 200/im in the rest-frame. Clearly this two-component 
model is not a unique description of the dust emission from 
SMMJ2135, but it does suggest that the system could be 
modelled by a combination of warm, compact and cooler, 
more extended dust component. 



3.1.2 



■ CO SLED 



To investigate the excitation within the ISM, we next con- 
struct the 12 CO SLED by calculating the 12 CO flux of each 
transition over the velocity range defined by the composite 
spectrum; we show the resulting SLED in Fig. [2] In Ta- 
ble 4 we also give the velocity-averaged brightness tempera- 
ture ratios between 12 CO(J up per) and 12 CO (1- 0), which we 
define as Rj.i = L'(J + 1 -> J)/L'(l 0) (|Greve et all 
l2003l ; iHarris et alj|2010h . The 12 CO SLED shows a contin- 
uous rise out to J up pcr = 6, followed by a sharp decline 



£ls ffir els Jnp 



10. The brightness ratios of all of the 



Jupper > 2 line luminosities compared to the 12 CO(1-0) are 
< 1, which interpreted in the framework of a single phase 
ISM would indicate that t he gas is sub-thermalis ed (the ther- 
malised prediction is ~ 1; iDevereux et al.lll994r i. Indeed the 
L'i2 CO ( 3 _ 2 j/L'i2 CO ( 1 _ ) line luminosity ratio for SMM J2135, 
0.68 ±0.03, is similar to that seen in local starb urst galaxies, 
0.64 ± 0.06 (IDevereux et al.lll994h . However, as IHarris et all 
(|2010h discuss and as we show later in H3.2I these ratios 
are more likely to indicate a multi-phase ISM as found in 
star-forming regions locally. 



The peak of the SLED at J u 



6 is similar to 



that s een in nearby starburst galaxies such as NGC 253 and 
M82 (|Bradford et alj|2003t IWeifi et al.ll2005l; iPanuzzo et all 
| 20ld) and AGN. e.g. M rk231 (|Papadopoulos et al l [20071 : 
Ivan der Werf et al.ll2010P ). However, the overall shape of the 
SLED suggests that SMM J2135 has proportionally stronger 
12 CO (1-0) emission (compared to the higher J U p pcr transi- 
tions) than any of these galaxies, indicating the presence 
of an additional low-excitation gas phase. Nevertheless, the 
bulk of the cooling, 60%, arises through the 12 CO Ju PP cr = 



5-7 line emission^ with a further 20% at 



'CO J UE 



> 



Ivan der Werf et al.l (|2010T i have recently constrained the 
12 CO in Mrk231 up to J U D per = 13 with Hersc h el, bu ilding 
upon the earlier study of IPapadopoulos et al.1 ((2007) , and 
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Figure 2. The integrated 12 CO SLED for SMMJ2135 showing that the SLED peaks around J up pcr = 6 (similar to M82, but with 
proportionally stronger 12 CO(l— 0)). The central panel shows the results of LVG modelling applied to the integrated SLED, which 
requires two temperature phases to yield an adequate fit. In this model they have characteristic temperatures of Tki n = 25 K and 60 K 
and densities of n = 10 27 and 10 3 ' 6 cm -3 respectively (we associate these two phases with a cool, extended disk and h otter, more 
compact, clumps). The right-hand panel shows a similar comparison but now using a PDR model l lMeiierink et al.ll2007l ; see 33.2. li 
which again requires a combination of both low- and high-density phases to adequately fit the integrated SLED. 



show that, unlike nearby starbursts (e.g. NGC 253 or M 82), 
the 12 CO line luminosity is roughly fiat across Ju PP er = 5- 
13. They calculate that only 4% of the total 12 CO line 
luminosity arises from the three lowest transitions, com- 
pared to ~ 20% in the three lowest transitions of our 12 CO 
SLED (and 43% in the Milky Way). The luminous high- 
Jupper emission they see is likely due to the presence of 
additional excitation phases in Mrk231, either high excita- 
tion PDRs or more likely a n X-ray dominated region (XDR) 
i|Spaans fc Meiierink! 2008). As we see later, we also find ev- 
idence for at least two phases of material in SMM J2135. 

In comparison to higher-redshift sources, we see that 
the 12 CO SLED for SMMJ2135 has a broadly similar 
shape to that seen in other high-reds hift SMGs, e.g. G N 20 
i|Carilli et al.ll2010T ) or SMM J16359 (|Weifi et al.ll2005l ). In- 
deed, recent progress with EVLA and GBT has started to 
provide 12 CO(1-0) detections of SMGs to complement th e 
earlier 12 CO(3-2) studies from PdBI (e.g. lGreve et al.ll2005h . 
Using these two transitions, we find that the i?3,i bright- 
ness temperature ratio we derive for SMMJ2135 (Table 4) 
is comparable to typical SMGs: 7?3, i = 0.55 ±0.05 (Ivi son et 
al. 2010c) and R 3 ,i = 0.68 ± 0.08 ^Harris et al.ll2010l ). sug- 
gesting that the CO SLED results we derive for SMM J2135 
may be applicable to the wider SMG population. 

To investigate the excitation of the gas reservoir within 
this galaxy in more detail we exploit the fact that the shape 
of the 12 CO SLED can provide information on the under- 
lying gas density and temperature distribution s and use a 
spher ical large velocity gradient model (LVG; Weifi et al.l 
120051 ) to attempt to fit the observed SLED. LVG techniques 
are the most widely used radiative transfer model which can 
account for photon transport when spectral lines are opti- 
cally thick, and can be used for efficiently solving the radia- 
tive transfer equatio n when the molec ule level populations 
are not thermaliscd (Bavct ct al. 2 0061 ). 

Thus we next use the LVG code to model the 12 CO 
SLED. The model assumes spherical symmetry and uniform 
kinetic temperature and density, a CMB background tem- 
perature of Tcmb = 2.73K(1 + z) which at z = 2.3 is ~ 9 K, 



the collision rates from iFlower] (|200ll ) with an ortho/para 
H2 ratio of 3 and a ratio of the CO abundance to the ve- 
locity gradient of CO jdvjdr = 10" 5 pc (kms -1 )" 1 . To fit 
absolute line intensities the model further uses the source 
solid angle size, Og, which can be expressed in terms of the 
equivalent radius of a face-on disk, ro = D Ay/fas /it). The 
model returns the temperature, density and size of the emis- 
sion region. 

Modelling the 12 CO SLED in this way, we find that a 
single temperature and density phase is unable to fit the 
SLED, and that two or more phases are required. This is 
consistent with the implication of the model of the dust 
SED and the sub-millimetre morphology: we need at least 
two temperature phases to adequately describe this sys- 
tem. This need for multiple phases to fit t he ISM in a 
high- redshift galaxy is unusual (however cf. ICarilli et al.l 
2010), although it is not unexpected given that such multi- 
phase ISM s are c ommonly requi r ed in local galaxies (e.g . 



Wild et al 



Mao et al.l 



19921: iGuesten et al.l [l993l: lAalto et all Il995l : 



2000; Ward ct al. 2003h . Thus, motivated by the 



fit to the dust SED we identify a two-phase fit to the 
SLED comprising four clumps and an extended phase, which 
provides an adequate fit (Fig. |2J. In this model, the low- 



excitation phase peaks at J up 



3 and is diffuse, with 



10 cm 3 and Tki n ~ 25 K. The denser, mildly-excited 



phase peaks around J u 



has T kill 



60 K and 



n ~ 10 cm and contributes ~ 60% of the total lumi- 
nosity over all the 12 CO lines (Fig. The low-excitation 
phase has a 12 CO SLED which is very similar to the in- 
ner disk of the Milky Way, while the warmer phase is well 
matched to the inner few hundred parsecs of the starburst 
region in NGC 253. Within this two phase model, the total 
gas mass is M 9as - (4.0 ± 0.1) x 10 10 Mq, a factor of ~ 3x 
the mass we determine from 12 CO(1-0), assuming a = 0.8, 
corresponding to an effective a = 2.0 (simlar to t hat de- 
rived for the nucleus of Arp 220; jScoville et al.lll997l )). From 
here on we adopt this value as the gas mass of our system, 
but note some of the differences if we had used the 12 CO 
line luminosity and a — 0.8, as is commonly assumed in 
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the analysis of high-redshift galaxies. The combination of 
this higher value for a and our estimate of i?3,i would in- 
crease the gas mass by a factor of ~ 4x over that estimated 
by assuming a = 0.8 and R^.i = 1 as usually assumed a t 
high redshift (see also llvison e t al. 2010i; [Harris et aTl boiO'). 
Adopting a ~ 2 results in a gas mass fraction in the cen- 
tral lkpc radius region of S MM J2135 of 75% (and closer 
to the 100% for typical SMGs. iGreve et al.ll2005h subject to 
uncertainties on the dynamical mass due to the unknown 
configuration of the gas within the system. 

3.1.3 Atomic Carbon 

We have also obtained strong detections of the [Cl]( 3 Pi — s> 
3 Po) and [Cl]( 3 P2 — > 3 Pi) emission lines, and from these 
we can begin to investigate the properties of the ISM. We 
start by noting the similar detailed morphologies of the 
[Cl]( 3 Pi -»■ 3 P ) and the 12 CO (3-2) and 12 CO (1-0) lines 
in Fig. [T] which suggests that the emission is arising from 
the same mix of phases with characteristic temperatures of 
T £ 40 K. Indeed observations of neutral carbon transitions 
have shown that the ratio of Lj CI](3p2 ^ 3Pi) / Lj OI]( a Pi _ > 3Po) 
provides a sensitive probe of the temperature of the in- 
terstellar medium at moderate densities. Both lines have 
modest critical densities (n cr it ~ 0.3-1.1 x 10 3 cm -3 ) and 
are therefore often thermalised in molecular clouds with 
n 10 3 cm -3 . The lines arise from states with energy lev- 
els Ti = 23.6 K and T2 = 62.5 K above the ground state, 
and thus their ratio is sensitive to the gas temperature if 
T gas <: 100 K. We derive L| CI](3p2 ^ 3Pl) / L'r CI]( 3 Pl ^ 3Po) = 
0.37±0.02, which is similar to that measured in nearby star- 
bursts and galactic nuclei such as M 82 and NGC 253 where 
L[nT](3p 2 ^3p 1 )/ L [ nT](3p 1 ^3p > ~ 0-33 (uBennett et ai] |l994l; 
IWhite et alii 19941 ; llsrael et al.1l 19951 ; iBavet et alj|2004h . but 
larger than typically found in the cooler, dense cor es of gi- 
ant molecular clouds (e.g. IZmuidzinas et all 119881 ) . To be 
compatible with previous studies we define R[ci] as the ra- 
tio of the [Cl]( 3 P 2 -> 3 Pi) to [Cl]( 3 Pi -> 3 P ) temperature 
integrated line intensities in K km s _1 . In SMM J2135 we de- 
rive an integrated line ratio of R[ci] = 0.88 ± 0.01. We can 
use this ratio to estimate th e excitation temperature T ox 
following Istiitzki et all (| 19971 ). T cx = 38.8 K/ln(2.11/P [CI] ), 
and derive T ex = 44 3 ± 1.0 K for SMM J2135. 

We also follow IWeifi et alj (|2005f > to derive the total 
mass of neutral carbon via: 

Mci = 1.902 x 10- 4 Q(T cx ) e Tl/T -L; ci](3Pl ^3 Po) (1) 

where Q(T CX ) = l+3e" Tl/T<!X +5e" T2/Tcx is the [Cl] partition 
function and T cx is defined above. Using T ox = 44.3 ± 1.0 K 
and our measured luminosity of the upper fine structure line, 
[Cl]( 3 Pi — > 3 Po), we estimate a carbon mass of Mrcjn = 
(9.1 ±0.2) x 10 6 Mq. Combining this with our gas mass esti- 
mated from the LVG analysis we derive a [Cl] abundance of 
M([Ci])/6M(H 2 )= (3.8 ±0.1) x 10~ 5 which is higher tha n 
the Galactic value of 2.2 x 10 -5 from iFrerking et all l|l989f) . 
Similarly, the ratio of the [Cl] and 12 CO brightness tempera- 
tures, Lj CI ] (3Pi ^ 3Po) /Li 2co(1 _ ) = 0.41±0.02, falls towards 
the upper end of the range 0.2±0.2 observe d in local galaxies 
l|Gerin fc Phillips! |2000| ; IBavet et al.ll2006h . A higher abun- 
dance of [Ci] rela tive to 12 CO is ex pected in regions with 
lower metallicity l|Stark et al.lll997l ) where 12 CO is photo- 



dissociated, in regions with a high ionisation fraction which 
may drive the chemistry to equilibrium at a high [Cl]/ 12 CO 
ratio, and in starburst sys t ems with a high co smic ray flux 
jPapadopoulos et al J 12004 : iFlower et ai1ll994h . 

We can also attempt to model the [Cl] emission within 
the framework of our two-phase LVG model. We find that 
we are only able to simultaneously fit the [Cl] luminosities 
and ratios if the carbon abundance is allowed to vary be- 
tween the two phases. To fit the observed R[ci] ratio and the 
[Cl] to 12 CO brightness temperatures the LVG model sug- 
gests [Cl/ 12 CO]~ 0.2 for the warm, dense phase, which is 
close to the Galactic value, [Cl]/ 12 CO ~ 0.13 (|Frerking et all 
1989), while for the extended, low-density phase we require 
an abundance of [Ci]/ 12 CO~ 3.0, which suggests a deficit 
of oxygen, and hence possibly a very low metallicity. Us- 
ing the correlat ion between [Cl]/ 12 C O and metallicity (12 
+ log[0/H]) in iBolatto et al.1 (|200Cj ). for the cool compo- 
nent we derive 12±log(0/H)=7.5 and hence a metallicity 
of ~Z/Z0=l/25. This is similar to that of the Sextans 
dwarf galax y, to the lowest m etallicity, young starbursts 
seen locally (|Brown et al.ll2008h . an d to the lowest metallic - 
ity galaxies found at high redshift (|Yuan fe Kewlevl f2009). 
However, the high [Ci]/ 12 CO rat io could also be due to a 
high cosmic ray flux in the galaxy |Papadopoulos et al.1 |2004| ; 
llsrael fc Baasll2OO20 . Indeed a factor of ~ 10 x enhancement 
of cosmic ray flux will yield the high [Cl/ 12 CO] ratio we ob- 
serve and this may not be unlikely given the high star forma- 
tion r ate density of the star-forming regions (|Swinbank et al.1 
2010). Moreover, an enhanced [Cl/ 12 CO] rati o has also been 
observed in the starburst galaxy NGC 253 (|Harrison et al.1 
ll995l ; lBradford et alJliool ) attributed to a high cosmic ray 
flux. 

Finally, the [Cl] detection s can be combined with the 
observations of [Cn]157.7 from llvison et all (|2010t ) to inves- 
tigate the ISM cooling. We can compare the cooling in each 
of our lines relative to the bolometric (8-1000/im) luminos- 
ity, to assess their importance. The fraction of the bolomet- 
ric luminosity in the [Cll] line is 0.24%, the rotational 12 CO 
lines ( 12 CO(1-0) to 12 CO(10-9)) contribute 0.09% and the 
[Cl] lines result in just 0.03%. Therefore, [Cn] dominates 
the cooling with the total cooling due to 12 CO and [Cl] rep- 
resenting ~ 50% of the cooling due to [Cll]. [Cn] is one 
of the brightest emission lines in galaxies and can account 
for 0.1-1% of the far-infrared lumino sity of the nuclear re- 
gions of galaxies (|Stacev et al.1 ll99lT ) ; our source lies well 
within that range. This situation is in contrast to compa- 
rably luminous systems in the local Universe, e.g. Arp220, 
where the [Cn] emission is just 1.3 x 10~ 4 of Lfir (Gerin 
& Phillips 1998). This most likely results from saturation 
of [Cll] in very high-density star-forming regions (Luhman 
et al. 1998). The proportionally stronger [Cll] emission we 
see in SMM J2135 may then reflect slightly lower densities 
and more extended star formation (Ivison et al. 2010c), or 
be due to the lower metallicity of the gas compa red to local 
ULIRGs jlsrael et al.lll996l ; lMaiolino et al.ll2009h . 

Overall, our observations of atomic carbon suggest that 
it arises from a cool phase within the galaxy and that inte- 
grated over the whole galaxy the [Cl]/ 12 CO ratio is higher 
than the Galactic value. However, we also find evidence for 
variation in the ratio within the system, with the cool, low- 
density phase having significantly enhanced [Cl], compared 
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Figure 3. The relation between far-infrared luminosity and 
HCN luminosity for galaxies. We compare our observations of 
SMMJ2135 with the correlation seen locally and find that this 
high-reds hift galaxy f ollows the tight correlation seen in local 
galaxies ( Juneau et al. 2009). In the main panel we plot obser- 
vations from a sample of 34 local galaxies (star- forming galaxies, 
star-forming g alaxie s with AGN and strong AGN) taken from 
lJuneau et al.l ll2009h . In the inset panel we compare our upper 
limit on the ratio of L^CN d— 0) / L'i2QQfi_p i to the ratio in local 
LIRGs and ULIRGs from iGao fc Solomonl [|2004h showing that 
this upper limit is consistent with the lGao* fe Solomonl ^2004) re- 
lation. Again we see that the properties of the dense gas and 
far-infrared emission in this galaxy are similar to those seen lo- 
cally. 

to 12 CO, suggesting lower metallicity in this phase or a 
higher cosmic ray flux in the overall system. 

3.1.4 HCN (3-2) 

Our spectra also cover emission lines from a number of 
species other than 12 CO; notably, we detect HCN(3-2) emis- 
sion (Table 2). HCN is an effective tracer of dense gas due 
to its high dipole moment, requiring ~ 100 x higher densi- 
ties for collisional excitation than 12 CO(1-0). Indeed, HCN 
is one of the most abundant molecules at densities n <; 
3 x 10 4 cm -3 (compared to critical densities of Z 500 cm -3 
for low-Jupper levels of 12 CO). We note that the velocity cen- 
troid of the HCN emission is redshifted by approximately 
+230 ± 100 kins -1 , relative to the nominal systemic redshift 
of the system derived from 12 CO. As we discuss below, this 
hints that the HCN emission may arise predominantly from 
only one of the kinematic components within the galaxy. 

It has been shown that the HCN line luminosity is 
tightly correlated with far-infrared luminosity in local spi- 
rals and ULIRGs, with a linear relation holding over at 
least three decades in luminosity, and a similar relation 
may hold at hi gh redshift (Fig. IGao fc Solomonl |2004 
IGao et all 120071 1. This suggests that HCN is a good tracer 
of the dense gas which fuels the massive star formation, 
which in turn is responsible for the far-infrared emission. 
In Fig. [3] we show the relation between far-infrared lu- 
minosity and HCN(3 -2) line luminosity for local spirals, 
LIRG s and ULIRGs l|Gao fc Solomonl |2004 I Juneau et ail 
120091 ). SMMJ2135 lies at the high- luminosity end of this 



correlation, corresponding to the highest-luminosity local 
ULIRGs, and is consistent with the relation. In the inset 
panel on Fig. [5] we also compare our s ource to local LIRGs 
and ULIRGs in lGao fc Solomonl |2004h . We show the upper 
limit on Lhqn(i-o) / L'i 2c0 (i_o) using our 3-er upper limit 
for HCN (1-0) flux, demonstrating that this is also consistent 
with the increasing Lh CN /L'i 2co ratio seen with increasing 
Lfir at low redshift. 

It is also pos sible to estimate a total gas mass from the 
HCN luminosity (|Gao et al .120051 ). Using our 3-cr upper limit 
on the HCN(l-O) luminosity we derive a limit on the dense 
gas mass of M dGnse (H 2 ) = a H CNri-o) Litr w < 4.5 x 10 10 M 
with ohcn < 10 l|Gao fc Solomonl 120041 ). We caution that 
there is at least a factor 3 x uncertainty arising from the un- 
certainty in ohcn- The value of ohcn we adopt is appropri- 
ate for a virialised cloud core; however, ULIRGs and LIRGs 
usually have a much higher HCN brightn ess temperature, 
resul ting in their «hcn being much lower (|Gao fc Solomonl 
12004) . and hence this is reported as an upper limit on the 
gas mass. Since this gas mass is consistent with the mass 
derived from our LVG analysis of 12 CO, this suggests that 
with better estimates of Qhcn, HCN may be a promising 
route to constraining the dense gas masses of high-redshift 
galaxies in the future. However, we note that the consistency 
of 12 CO-based and HCN-based gas mass estimates may be 
coincidental if these lines trace different phases. We further 
caution that at least two HCN lines may be required to de- 
rive the correct value of qhcn due to the strong variations in 
HCN(4-3) /H CN(l-O) and HCN(3- 2) /HCN(l-O) luminosity 
ratios locally (|Papadopoulosi 120071 ) . 

As with [Cl], we also attempted to use our two-phase 
LVG model from the 12 CO SLED to model the HCN(3-2) 
and HCN(l-O) emission from this system. However, we find 
that any models which fit the HCN SLED, in particular the 
bright HCN(3-2) emission, over-predict the luminosities of 
the 12 CO SLED and result in gas masses which are greater 
than the dynamical mass. A similar problem has been iden- 
tified in local ULIRGs which exhibit enhanced HCN line 
emission and high ratios of Lh CN /L'i 2co due to a higher 
proportion of dense gas in the s e strong starbur st systems 
dSolomon fc Vanden Boutll2005l ; IGao et aJj[2007h . 

3.1.5 H 2 

Finally we comment on the limits obtained for H2O. H2O 
emission arises from the warm, dense gas in the densest 
regions of a starburst or around an AGN. The strength 
of the emission is thus a probe of the radiation density 
and so reflects the compactness of the far-infrared source. 
The H 2 (21,1-20,3) (731.681GHz) emission line has re- 
cently been detected in the local ULIR Gs Mrk231 and 
Arp220 (e.g. iGonzalez- Alfonso et alj|2010h revealing warm, 
dense material, possibly ari sing from an XDR associ ated 
with AGN activity (see also ISpaans fc Mei ierink 200^ and 
van der Werf et all [2010). We compare our limits on the 
H2O line luminosities (as a fraction of bolometric luminos- 
ity) with Mrk231 and Arp220, where the luminosity ra- 
tios are Lh 2 o/Lpir ~ 1-2 x 10~ 6 . We derive limits on the 
line luminosity ratio of Lh 2 o/Lfir < 1-5 x 10 -6 from the 
H2O (5i,5-4 2 ,2) and H2O (21,1-20,2) transitions. These sug- 
gest that if an XDR exists within SMM J2135 it is unlikely 
to be more luminous than those seen in these local AGN- 
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dominated ULIRGs. This conclusion is reinforced by the 
relatively low luminosity of the high-J U p pcr 12 CO lines in 
SMMJ2135 compared to the substantially brighter high- 
J lines in Arp220 and Mrk23 1 which again require an 
XDR to fit their 12 CO SLEDs (|Spaans fc Meiierinkl [200I : 
Ivan der Werf et al.ll201ot ). 

3.2 Integrated Properties: Physical Properties of 
the ISM 

When massive stars are formed around a molecular cloud, 
their ultra-violet (UV) radiation changes the chemical prop- 
erties of the cloud's surface layers; the far-UV flux pho- 
todissociates the outer layer and ionises the material (such 
as carbon). This layer then cools primarily via atomic fine 
structure lines of [Ol], [Cn], [Ci], and the rotational CO 
lines. Emission from these differe nt coolants arises fro m dif- 
ferent depths within the PDR (|Kramer et al.l [20041 . such 
that the surface layers are dominated by emission from hy- 
drogen, [C11] and oxygen, although further into the star- 
forming clouds (as the extinction, Ay, increases), hydrogen 
becomes molecular, whilst the ionisation of carbon declines, 
with [Cn] becoming [Ci] and then combining into CO. More- 
over, with increasing density, the higher- J uppcr 12 CO emis- 
sion becomes stronger, as do other molecular ga s tracers 
such as HCN, CS and CN jKaufman et al.lll999l ). Deeper 
into the cloud, the gas is molecular, but still has a higher 
temperature than in the far-UV shielded core. Significant 
effort has been devoted to modelling (and predicting) the 
fine structure and molecular emission line ratios associated 
with these star-forming regions using PDR models. These 
account for variations in density, te mperature, dumpiness 
and time-dependent chemistry (e.g. iMeiierink et al.l 120071 : 
iKaufman etai1ll999f ) and using their predictions of the ra- 
tios of [Cn], [Ci] and 12 CO we can investigate the typical 
far-UV intensity and ISM density within the star-forming 
regions in SMM J2135. 

3.2.1 PDR models 

To model the full-range of observed li ne ratios from 
SMM J213 5 we use the PDR m odels of IKaufman et all 
l | 19991 ) and IMeiierink et all (|2007h in which the emission is 
determined by the atomic gas density (n; the density of H 
nuclei) and the incident far-UV radiation field from mas- 
sive stars ijnv — 6-13.6 eV) which is expressed in terms 
of Go (where Go is the average far-UV radiation field in 
the Milky Way and is expressed in Habing units: 1.6 x 
10 -3 ergs -1 cm -2 ). Simple PD R models (including one of 
the two models we compare to; IKaufman et al . 1999) calcu- 
late the line emission generated by a cloud illuminated only 
on one side. However, for observations of external galaxies 
it is not clear that such a "single cloud" model will ade- 
quately describe the observed molecular lines. For clouds 
illuminated on all sides we note that an observer will de- 
tect optically thin radiation emitted from both the near and 
far sides of the cloud, but only see optically thick emission 
from the near side. Generally, low-J uppe r 12 CO transitions 
are assumed to be optically thick f or the nominal Ay = 10 
cloud depths in the PDR models (jHailev-Dunsheath et al.l 
l2010bl ). whilst [Cn], [Ci], 13 CO and L F ir are assumed to 



be optically thin (although note our earlier estimates of the 
high column densities and extinction in this system). We 
have not applied any corrections to the emission lines to 
account for these differences between optically thick and 
optically thin emission. Instead, we also compare to more 
complicated models involving two or more clouds, with dif- 
fering densities and incident rad i ation fields which have been 
constructed by IMeiierink et all (2007). The implementation 
of these models assumes all emission lines are coming from 
the same gas in the same region, which we know to be a 
vast over-simplification; however, these solutions do provide 
us with order of magnitude estimates of the characteristic 
conditions in the system. 

In Fig. [4] we show the model grid of n and Go val- 
ues for various emission line ratios for SMMJ2135. The 
12 CO J/(J - 1) and [Cl]/ 12 CO line luminosity ratios vary 
over six orders of magnitude in Go; however, the [Cii]/FIR 
and [Cll]/ 12 CO ratios appear to break this degeneracy. 
Indeed, in the n ~ 10 3 -10 5 cm -3 regime, these provide 
a strong constraint on Go. Using the two different PDR 
models we convert the allowed parameter space for each 
line ratio to derive the peak likelihood solution, which is 
shown in Fig. [4] for both models. We find that the solu- 
tions cluster around moderate densities (n ~ 10 4 cm -2 ), 
similar to those claim ed for local ULIRG s from PDR mod- 
elling of H2 emission (|Davies et all 120031. The peak likeli- 
hood solutions use all the line ratios available aside from 
the Lj CI]( 3p 2 _ + 3p l) /L{ CI](3 p i ^3p o) ratio. This line's track 
deviates strongly from the preferred solution from the other 
lines. We expect this reflects the sensitivity of this line ratio 
to details of the PDR, such as geometry, and we discuss spe- 
cific problems with the PDR model treatment of [Ci] below. 

The best-fit solution to our suite of line ratios is n = 
10 4 1±0 3 cm - 3 and G = iq 3 - 6 ±°- 7 H abing fields for the PDR 
models from IKaufman et all (|l999l ) and n ~ iq 4 - 3±0 - 2 C m -3 
and G o ~ 10 3 ' 3±0 2 Habing fields from IMeiierink et al.l 
(2007). The best-fit solutions are within 1-a and are in 
agreement with earlier estimates ba sed on just 12 CO (1-0), 
[Cn] and Lfir in llvison et al. (2010). As we see from Fig. [H 
these physical co nditions are similar to local starburst g alax- 
ies and ULIRGs ijStacev et alJll99ll : lDavies et al]l2003 ) and 
at the upper end of both density and Go for normal star- 
forming galaxies ()Malhotra et al.ll200ll ). for example the ra- 
diation field is ~ 1000 x more intense than in the disk of 
the Milky Way. The characteristic density compares well to 
the dense star-forming cores of Galactic molecular clouds, 
or to the central compact nuclear gas disks of local ULIRGs 
such as Arp220 |Downes fc Solomon! [T998; Saka moto et al.l 
[l999j) and is also consistent with the integrated 12 CO SLED 
(see £|3.1.2[) . This characteristic density we derive is greater 
than the critical density for [Cn], [Ci] and CO(1-0), above 
which the intensity of the lines saturates. Similarly, the Go 
value we derive suggests the condition s in the PDRs are a t 
the limit of the steady state solution l|Kaufman et alj|l999t ): 
with radiation pressure on grains comparable to turbulence, 
suggesting that the PDRs may be ionising material. 

Our results imply that the underlying physics of the 
star formation occurring in the system is similar to the dense 
star-forming cores of Galactic molecular clouds (character- 
istic of massive OB star- forming regions), even though the 
total energetics of the system is f ar more extrem e than in the 
disk of our galaxy. Interestingly, iTevssier et al.l (|2010T ) have 
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Figure 4. Left; Luminosity line ratios (in units of Lq) from 12 CO, [Ci] and [Cn] as a function of density and far-UV flux (Go in units of 
Habing field), from the PDR models of iKaufman et all il999T ). Tracks are drawn for the measured ratios in solid lines with the 1-a errors 
as dotted lines. The line ratios within SMMJ2135 inte rcept at n ~ 10 41 cm" 3 an d Go ~ 10 3 ' 6 Habing fields and n ~ 10 4 ' 3 cm -3 and 
Go ~ 10 3 ' 3 Habing fields for the Kaufman et al. | \\m$) and lMeiieriiik" et al. models respectively. To better display the preferred 

regions of parameter space we combine the probability distributions from all of these lines (excluding £[ci](2— i)/-^[Cl](l— 0) f° r reasons 
explained in § 13.211 to derive a peak likelihood solution for each model. Right: The same parameter space for the models, now comparing 
the peak likelihood solutions from the model grids to the derived values of Go and density (n) for various low- and high-redshift galaxies 
and molecular clouds. We also indicate the regions of parameter space which ar e typically encompa ssed by Galactic OB star-formation 
regions, starburst galaxies and non-starburst and Galactic molecular clouds from lStacev et aU l ll99ll ). This shows that the line ratios for 
SMM J2135 are consistent with those typically found in local starbursts or ULIR Gs. The local galaxy that lies within the contours of the 
peak likelihood for the Kaufman models is NGC 253, a sp i ral st arburst galaxy l|Negishi et al]|200lh . We also show the peak likelihood 
solution derived from the PDR models of Mciicrin k et alj d2007f) which provide a concomitant solution. The contours represent the 1- 
and 2-<t limits compared to the best-fit solution. The peak likelihood solutions are derived without including the L[ci](2— 1) /L[ci](i-o) 
track, which we plot on both panels to show its discrepancy from the general solution from the other lines (see § 13. 2D . 



recently used hydrodynamic simulations of major mergers 
of disk galaxies to show that the density of gas at which 
a majority of the star formation in these systems occurs is 
n~ 10 4 cm -3 , as we find here. 

Of all the line ratios we analyse, the most obvi- 
ous outlier in the PDR modelling is the [Cl]( 3 P2 — >■ 
3 Pi)/[Cl]( 3 Pi ->■ 3 P ) emission line ratio (Fig. Q}. How- 
ever, we note that this line ratio is one of the most sensitive 
to Go and n and it would only require a modest increase 
in the predicted [Cl]( 3 Pi — > 3 Po) flux in the models, ~ 3x 
stronger, to shift this track into agreement with the solution 
based on the bulk of the other ratios at n ~ 10 4 cm -3 and 
Go ~ 10 3 ' 5 . This problem may arise because the PDR mod- 
els assume a particular geometry, which determines the ratio 
of diffuse to dense gas, and assume a homogeneous medium 
(an unrealistic assumption for our source). The resulting lu- 
minosity ratios based on 12 CO/[Cl] and [Ci] are sensitive to 



this geometry and the [Cl]( Pa 



3 P 1 )/[Ci]( 3 P 1 -> 3 P ) 



ratio is particularly s ensitive to the homogeneity of the 
medium (jSpaansI I1996I K as demonstrated by the difficulty 
the models have in reproducing the [Cl]( 3 Pi — ¥ 3 Po) or 
[Cl]( 3 p2 — ¥ 3 Pi) emission lines even in tranquil environ- 
ment s (e.g. iPineda fe Bensch 2007; iPapadopoulos fe Grevd 
l2004l h Indeed lPapadopoulos fc Grevd (|2004T ) argue that [Ci] 
is distributed relatively ubiquitously in molecular clouds 
(the critical density for excitation of both [Cl]( 3 Pi — s> 
3 Po) and [ClK 3 P 2 ->■ 3 Pi) is roughly that of 12 CO(1-0); 
ISolomon fc Vanden Boutl [2005 V This is contrary to PDR 
theory, which places [Ci] in a narrow [Cn]/[Cl]/ 12 CO tran- 



sition zone on the surface of far-UV illuminated molecular 
clouds, and potentially makes [Ci] a very effective tracer of 
H2 mass. Therefore our low ratio may be arising from ge- 
ometrical optical depth effects in which the [Ci] effectively 
follows the volume (rather than the PDR surface), which 



would naturally lo wer the [ClK P 2 



3 Pi) /[Cl]( 3 P! 



3 Po) ratio (e.g. see iKramer et all l200fet ) as would happen 
from any process which enhances [Ci] emission through the 
cloud's volume (e.g. cosmic rays). We also caution that the 
[Cl]( 3 P 2 3 Pi) /[Cl]( 3 Pi -»■ 3 P ) ratio has been m easured 
in only a few local galaxies (e.g. Bayet et al j 2006 ) and so 
the fine structure excitation of neutral carbon has yet to be 
fully characterised and calibrated locally. 



3.2.2 Limits on 13 CO, CN, HNC and HCO+ 

Our observations have also simultaneously covered a number 
of dense gas tracers and from these we can place limits on 
the line luminosities of 13 CO, CN, HNC and HCO+. We 
can then ask if these limits are consistent with our preferred 
characteristic density and UV radiation field from the PDR 
model. 

First we consider the limits on the 13 CO. The 13 C nu- 
clei are synthesised by CNO processing of 12 C nuclei from 
earlier stellar generations ( 12 C is produced from He burn- 
ing on rapid timescales in massive stars). Since 13 C is pro- 
duced more slowly in intermediate mass stars during the 
red giant phase (at ~ lGyr), the 12 C/ 13 C ratio has been 



© 0000 RAS, MNRAS 000, 000-000 



12 Danielson et al. 



Table 3. 
Model Kinematic Parameters 



Component 


V 


a 




(km s _1 ) 


(km s" 1 ) 


z 2 


-167 ±9 


75 ±8 


Zx 


28 ±9 


75 ±8 


Y 


165 ±13 


157 ± 17 


X 


396 ±9 


76 ±9 



Notes: Velocities and line widths of the three kinematic compo- 
nents seen in the CO and [Cl] spectra. The velocities are given 
with respect to a heliocentric redshift of z = 2.32591. 

proposed as a diagnostic of the nucleosynthesis history in 
galaxies. A high ratio would then be seen at high redshift as 
the young starbursts will have had insufficient time to form 
secondary nuclei of 13 CO. In the Milky Way and nearby 
star-forming galaxies, the 12 CO/ 13 CO line luminosity ra- 
tio (L'i2 CO /L'i3 CO ) is ~ 5-10. However, local ULIRGs have 
shown a deficiency of 13 CO (possibly due to gas inflow; 
Rupke et al. 2008), which raises the 12 CO/ 13 CO line ratio 
to ~ 30 (Greve et al. 2009). Indeed, recently, Henkel et al. 
(2010) showed that the 12 CO/ 13 CO luminosity line ratio in 
the Cloverleaf quasar (0=2.6) is L/ 12co /L/i 3co =40±17, sug- 
gestive of a deficiency of 13 CO. Within SMM J2135, we esti- 
mate limits on the ratio of >5.8 and >6.5 for the J up por=l 
and J uppC r=3 transitions respectively. These limits are con- 
sistent with the Milky Way, but more sensitive observations 
are required to place firm conclusions as to whether a defi- 
ciency exists similar to that observed in other local ULIRGs 
and high-redshift galaxies. 

We can compare these limits to the predictions from 
the PDR models of Meijerink et al. (2007). The model pre- 
dictions are Li2 CO /Li3 C o = 2-13 and 2.08 for the J U pp C r=l 
and J up pcr=3 transitions respectively, indicating lower 13 CO 
abundance in SMM J2135-0102 than the model predicts. 
This may arise due to a number of effects. It is known that 
the line luminosity ratios of giant molecular clouds are lower 
than galaxy-integrated properties due to the lower opacity 
of 13 CO compared to 12 CO in the inter-cloud medium. Since 
we are measuring galaxy integrated properties, this would 
naturally lower the ratio. Clearly it would be useful to ob- 
tain a more sensitive observation of the 13 CO to determine 
the extent of this discrepancy with PDRs. 

Finally, we note that the line luminosity ratio of the 
dense gas tracers of CS(2-1)/ 12 CO(1-0)<0.13 is lower than 
that predicted from the PDR model (0.43). In contrast, 
the HNC(3-2)/ 12 CO(l-0)< 1.3 and HCO+(3-2)/ 12 CO(l- 
0)< 1.2 are consistent with those predicted by the PDR 
models (1.3 and 0.15 respectively). A more detailed test of 
the model prediction for these dense gas tracers will again 
require more sensitive observations. 



3.3 Kinematically Resolved Properties: 
Decomposed 12 CO SLED 

With the high signal-to-noise of our observations of 12 CO 
and [Cl] in Fig. [T] it is obvious that the spectra show sig- 
nificant kinematic structure and that more intriguingly this 
structure appears to vary between the different transitions 




Figure 5. The 12 CO SLEDs for the three kinematic components 
seen within the 12 CO spectra (from left-to-right for the compo- 
nents at X, Y and Z at ~ —180, and +200 km s -1 respectively). 
We overplot on each of these a scaled version of the best-fit LVG 
model from the combined SLED. This shows that two compo- 
nents, X and Z, also peak at around J up pcr ~ 6, while compo- 
nent Y appears to peak slightly lower, Juppor ~ 4, suggesting that 
there may also be temperature differences between the kinematic 
sub-components within the galaxy. Interestingly, component Y 
also appears to contain the bulk of the cold gas as traced by 
12 CO(1-0) compared to either X or Z. 
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Figure 6. The correlation between measured gas and SFR den- 
sities for SMM J2135 comp ared to samples o f local star-forming 
galaxies and ULIRGS from iKennicud jl998h and high-redshift, 
as derived from 12 CO(l— 0). The dashed line is the best fit power 
law to local data with E S F oc S^ 27 (Genzel et al. 2010). We 
see that the star formation within SMM J2135 is consistent with 
this relation assuming Tgas = r S p = 1 kpc. SMM J02399-0136 
is taken from Ivison et al. (2010), GN20 from lCarilli et al]|2010l 
and the BzKs from lAravena et al. which are unresolved 

in 12 CO(1-0), hence the arrows representing the lower limits in 
measured gas density. 



and species. We therefore attempt to fit a model to the 12 CO 
spectra and extract 12 CO SLEDs for each kinematic compo- 
nent, in order to search for excitation structure within the 
system. 

To construct a kinematic model for the 1 CO emission 
we fit the composite spectrum (Fig. [1} with a series of in- 
creasingly complex models. We find that a double Gaus- 
sian profile fit provides a ~ 14-rr improvement over a single 
Gaussian profile, with a triple profile giving a further ~ 6- 
cr improvement over the double Gaussian profile. Although 
this model reasonably describes the line profile of the sys- 
tem, P. Cox et al. (in prep) use the velocity structure from 
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Table 4: Line ratios 



Component 


Mu 2 +He 

(1O 9 M ) 


no, 
^3,1 




^5,1 




R 7 ,l 


Rs,i 


CI(2-1) 
r I 

CI(l-O) 


CI(l-O) 
jl 

CO(l-O) 


Total 


40.0 ±0.3 


0.68 ±0.03 


0.50 ±0.04 


0.35 ±0.02 


0.28 ±0.02 


0.119 ±0.008 


0.064 ±0.005 


0.37 ±0.02 


0.41 ± 0.02 


Z 


10.1 ± 1.8 


0.71 ±0.19 


0.42 ± 0.12 


0.51 ±0.11 


0.41 ±0.09 


0.17 ±0.04 


0.09 ±0.01 


0.85 ±0.21 


0.32 ± 0.10 


Y 


23.1 ±3.2 


0.67 ±0.14 


0.49 ± 0.10 


0.22 ±0.05 


0.19 ±0.05 


0.09 ±0.02 


0.04 ±0.01 


0.19 ±0.05 


0.43 ± 0.09 


X 


6.8 ± 3.2 


0.70 ± 0.48 


0.67 ±0.42 


0.48 ±0.29 


0.38 ±0.23 


0.06 ±0.06 


0.11 ±0.06 


0.43 ±0.27 


0.38 ± 0.28 



Masses and line luminosity ratios for individual kinematic components. 



*Rj\ represents L' 



CO(J+1-J)/ L CO(1-0) 



high-resolution interferometry with PdBI up to J UF 



7 



associate with the core of the second interacting system, is 



to derive a map of the dynamics of SMM J2135 on ~ 100-pc 
scales. The system can be described by two interacting disks: 
the lower-redshift of these displays a double-horned profile 
(termed Z\ and Z2), with a second, Y, at higher velocity 
which appears to connect to a further structure, X, extend- 
ing to much higher velocity. Based on this interpretation 
we construct a model which comprises three components: 
two single Gaussians representing X and Y and a double- 
Gaussian for Z (constrained to have equal line widths and 
intensities for the two components). We fit this model to 
the composite spectrum, allowing the intensity, width and 
velocity of each component to vary and show the resulting 
best fit to the average line profile in Fig. [T] and report its 
parameters in Table 3. This fit has a x 2 of ~ 20 for seven 
degrees of freedom and provides an adequate description of 
the composite line profile. 

Using the best-fit model parameters from the composite 
12 CO spectrum, we apply the fit to each 12 CO spectrum in 
turn, in order to extract the SLEDs for the three components 
in the model. Beginning with the average model parameters, 
we perturb the best fit randomly by up to ±20kms~ 1 and 
fit this new realisation to each of the 12 CO spectra in turn, 
allowing only the relative intensities of the components to 
vary. We then sum the \ 2 f° r the fit to each line for each 
model and determine the parameters for which the total \ 2 
is a minimum. We show the best fit model's line profile on 
each line in Fig. Q] and report their parameters in Table 3. 
To calculate errors, we allow the best fit parameters of the 
kinematic model to randomly vary (by their l-er uncertain- 
ties) and run 10 Monte Carlo simulations, determining the 
1-cr variance in the range of fluxes from those models where 
the total x 2 of the model fit is within Ax 2 = 21 of the best 
fit (since we have seven emission lines and fit with three 
free parameters to each line). The errors on the parameters 
for the model fit to the average (velocity and dispersion in 
Table 3) are calculated using the standard deviation of the 
parameters for all models for which the total \ 2 of the model 
fit is within Ax 2 = 21 of the best fit; these are quoted in 
Table 3. 

The 12 CO SLEDs for the three individual components 
are shown in Fig. [S] We derive the gas masses for each com- 
ponent using the 12 CO (1-0) emission line and assuming our 
effective conversion factor a — 2.0 and list these in Table 4. 
The most highly excited components are X and Z, peaking 
at Juppcr ~ 6, and containing ~ 15% and ~ 25% of total 
gas mass respectively (Table 4). Component Y, which we 



somewhat less excited, peaking at J u 



4, and contains 



most of the gas in the system (~ 60% of the total gas mass 
and a majority of the cold gas). 

For each of the SLEDs, we perform further LVG mod- 
elling as applied to the combined spectrum. However, within 
the errors, we cannot measure any temperature differ- 
ence between the three components although there is a 
hint that Y is cooler than X or Z. However, while the 
[Cl]( 3 P 2 -> 3 Pi) and [Cl]( 3 Pi ->■ 3 Po) lines have lower 
signal-to-noise than the 12 CO lines, we can also decompose 
them using the same method. Table 4 shows the ratio of 
L( C i](3p 2 ^ 3 Pl) / Lj CI]( 3 Pl _> 3 Po) for the different components, 
with Z having the largest ratio and Y having the smallest. 
The corresponding excitation temperatures are T cx > 100 K, 
T cx = 24.7±0.5K and T cx = 52.5±6.3K for Z, Y and ^re- 
spectively (where the temperature calculated for the Z com- 
ponent may be discrepant and likely caused by a problem 
with the model fit to the line; see Fig. 1). This supports the 
suggestion from the 12 CO SLED that both X, and especially 
Z, are hotter than the more massive Y. 



3.4 Physical Interpretation 

ISwinbank et~aH (|201Ch show that the rest-frame ~ 250-/im 
emission from SMM J2135 is concentrated in four bright, 
star-forming regions each of which has a radius of ~100pc 
spread across a more diffuse structure with a total extent of 
2kpc in projection. Following this we fitted the integrated 
dust SED of the source (see H3.1.ip assuming that ~ 20% of 
the 870- /^m flux is coming from a relatively cool (Td ~ 30 K) 
extended dust component with a predicted radius of ~ 1 kpc 
and ~ 80% emerges from the four ~ 100-pc radius, hotter 
(Td ~ 60 K) and optically thick, clumps. Our modelling of 
both the 12 CO SLED and the kinematically decomposed 
[Cl] line ratios supports the presence of two phases in the 
ISM within this galaxy with a cool, 15-25 K, phase and a 
hotter, 45-60 K phase (although it is likely that these ac- 
tually represent a continuum of properties). We also find 
an enhancement in the [Cl]/ 12 CO ratio (likely due to low 
metallicity and/or higher cosmic ray flux) which is asso- 
ciated with the cool phase in this system. Our kinematic 
model identifies a component Y which is relatively cool and 
contains the about half of the gas in the system, while the 
hotter material is more closely tied to component Z (and 
to a lesser extent X). It is tempting to associate these kine- 
matic components with the clumps seen in the SMA map, 
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but as we discuss below, it is more likely that the clumps 
lie within these kinematic structures, rather than having a 
unique one-to-one relationship. 

Our observations also provide other limits on the char- 
acteristic size and morphology for the emission regions. 
We can use the characteristic density derived from the 
PDR analysis, together with the total gas mass to derive 
an effective radius of the system. Adopting M(H2+He) ~ 
4 x 10 10 Mp) and a density of ~ fO 4 cm" 3 from the 
iKaufman etafl (|l999h models, we derive the radius of an 
equivalent sphere of ~ 300 pc (or ~ 700 pc for a disk with a 
thickness of ~ 100 pc). These sizes are larger than the clump 
sizes in the sub-millimetre, which again is consistent with 
the interpretation from the LVG modelling that the 12 CO 
is much more wide spread than the star-forming clumps. If 
we similarly split the gas mass into four equal clumps, then 
the expected sizes of these are ~ 200 pc in radius, which is 
slightly larger than the far-infrared sizes of the star-forming 
regions within SMMJ2135. 

We can also use the value of Go ~ 10 3 ' 6 Habing fields, as 
determined from the PDR models, to estimate the character- 
istic size of the galaxy. Assuming a Salpeter IMF and a con- 
stant SFR in the burst then the ionising flux (6-13. 6eV) is 
~ 50% of the bolometric luminosity for ages of 10-100 Myrs. 
We can determine the size of a region where the typical 
ionising field, as estimated from our bolometric luminosity, 
equals Go. This gives an estimate for the system's radius of 
~ 1 kpc, similar to the characteristic size inferred from the 
LVG modelling of the dust SED and the overall extent of 
the system from the SMA observations. 

We can use the kinematics of the components and the 
potential sizes and estimate the gas mass fractions in these 
various components. We first note that the dynamical mass 
of a uniform sphere with velocity dispersion of a ~ 75- 
160kms _1 (Table 3) and an effective radius of r ~ 200- 
lOOOpc is ~ 0.1-3 x 1O 1O M . The gas masses from Fig. [5] 
(using a — 2.0) are ~ 7-23 x 10 9 M©, similar to the range 
of dynamical masses, indicating that the gas reservoirs in 
these components must have effective radii comparable to 
the estimates above from the various size indicators. This 
also shows that we cannot significantly increase a, for ex- 
ample to the Milky Way value of a ~ 4.6, since this would 
result in a range of gas masses of ~ 1.4-4.7 x 10 10 Mq, in 
excess of the dynamical masses unless the gas has disk-like 
kinematics and we are seeing these almost face on, which is 
unlikely given the large velocity dispersion. 

Our size estimates for the gas reservoirs are larger than 
the size of the clumps in the SMA map and are more compa- 
rable to the extent of the whole system, which suggests that 
the cold gas is more widely distributed than the g 100-pc ra- 
dius far-infrared clumps. Thus, these clumps are likely to be 
knots of emission lying within a colder and more extended 
structure. Although the clumps will trace the kinematics of 
the gas in which they are embedded, they may not represent 
the deepest parts of the potential well, and so the star forma- 
tion may be free to migrate around the gas reservoir. These 
intense star-forming regions have crossin g times of just a 
few M yrs and they have free-fall SFRs (|Krumholz fc Tar] 
120071 ) of ~ 0.001 indicating that they are supported against 
collapse (most likely by internal turbulence). 

Finally, with an estimate of the overall extent of the 
gas reservoir, we investigate where SMMJ2135 lies on the 



Kennicutt- Schmidt relation, which links the star-formation- 
and gas surface density (Esfr and E gas respectively). This 
relation is approximately linear up to a threshold of ~ 
O.IMqpc -2 , above which the density of the ISM becomes 
dominated by molecular ga s and the star-fo rmation law in- 
stead follows Esfr oc Eg a t (|Kennicutt|[l99 8) . Previous work 
has shown broad universality of this relation with redshift, 
although at z ~ 2, both surface brightness and instrument 
limitations mean that a direct comparison to low-redshift 
galaxies has been difficult, with extrapolations from high- J 
CO in intri nsically luminous g alaxies providing most con- 
straints fe.g. lGenzel et af1l2010h . Some of this tension can be 
alleviated by using the CO (1-0) emission (or LVG modelling 
of the full CO SLED) to provide a more reliable estimate 
of the gas mass and spatial extent, and in Fig. [6] we show 
the Kennicutt-Schmidt relation and overlay the position of 
SMMJ2135. We also include on the plot the high-redshift 
galaxies where CO (1-0) observations are available. This fig- 
ure shows that SMMJ2135 lies at the high star-formation- 
and gas-density end of the correlation, although within the 
scatter seen for local ULIRGs. 

Although we do not have spatially resolved spec- 
troscopy to probe the gas properties on ~ 100 pc scales, we 
can ask where the individual star-forming regions would lie if 
they too followed the Kennicutt-Schmidt relation. Assuming 
that 80% of the bolometric luminosity arises in these clumps, 
and that each clump has a radius of rsF=r gas ~ 100 pc, 
then the inferred gas density would be n ~ 3 x 10 5 cm -2 , 
which is substantially higher than that inferred from the 
PDR and LVG modelling of the dense gas tracers (n ~ 
i()3.5-4.5 cm -2 ). However, taking the remaining 20% of the 
bolometric luminosity with a radius of 1 kpc (and assuming 
the gas lies in an extended disk), the density is n ~ 0.2- 
1 x 10 4 cm 2 (with the main uncertainty being the disk thick- 
ness). Thus, although this calculation should be considered 
crude, the results are consistent with a model in which 
the sub-millimetre emission predominantly traces the star- 
forming clumps on ~ 100 pc scales, but these clumps are 
simply embedded in a much more extended structure which 
dominates the 12 CO emission. These results may also imply 
a break-down of the Kennicutt-Sch midt relation on scale s of 
individual GMCs (< 100 pc, see also lOnodera et alj2oToh al- 
though clearly high resolution millimetre observations would 
accurately constrain their sizes, which is a fundamental step 
in understanding the conditions within the ISM of high- 
redshift galaxies (P. Cox et al. in prep). 



4 CONCLUSIONS 

We have presented detections of a number of rest-frame, 
far-infrared and sub-millimetre molecular and atomic 
emission lines in the spectrum of the lensed sub-millimetre 
galaxy SMM J2135. We have used these lines to construct 
the 12 CO SLED up to Juppcr = 10 and analysed this and 
the other species, including the HCN(3-2) and [Cl] fine 
structure lines, to investigate the physical conditions within 
the interstellar medium in this galaxy. Our main results 
are: 

The 12 CO SLED resembles that of local starburst galax- 
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ies. An LVG analysis of the 12 CO SLED shows at least two 
phases with temperatures T ~ 25 and 60 K. 

The dust SED is best fit by a model in which the galaxy 
comprises an extended component with radius ~ 1-2 kpc 
with Td ~ 30 K and four high-density, star-forming regions 
with sizes of order r ~ 200 pc and T<j ~ 60 K. These charac- 
teristics are roughly consistent with the rest-frame ~ 250- 
/im morphology of the system as seen by SMA and the two 
components may link to the phases identified from the 12 CO 
SLED. At Juppcr = 1-0, the cooler, diffuse phase appears to 
comprise 50% of the molecular gas. We also derive a car- 
bon mass of Mci = (9.1 ± 0.2) x 10 6 M and hence a carbon 
abundance of ~ 4 x 10 -5 which is slightly above the Galactic 
value. LVG modelling suggests that the dense, star-forming 
regions have abundances similar to that seen in the Milky 
Way, whilst the extended diffuse gas has lower metallicity 
(or a higher cosmic ray flux). 

Using the 12 CO emission and assuming a — 0.8 we esti- 
mate the total molecular gas mass from the 12 CO(1-0) emis- 
sion as M(H 2 ) = (1.4±0.1) x 10 10 M . Comparing this with 
that derived from Juppcr > 3 transitions, assuming ther- 
malised gas, shows that these predicted gas masses are a 
factor 1.5-3.6X too low, indicative of multiple gas phases, 
as seen in local starbursts. However, our LVG analysis in- 
dicates a higher gas mass, M(H 2 )~ 4 x 1O 1O M , yielding 
an effective CO-to-H2 conversion factor of a ~ 2.0. If these 
conditions are relevant in the wider SMC population then 
this suggests that previous, J up pcr > 3 12 CO observations of 
SMGs may have underestimated the molecular gas content 
of SMGs by at least a fa ctor of ~ 2-5 x (see also llvison et al.l 
l20ld ; lHarris et al.ll2010T l. 

Within SMMJ2135 we detect the HCN(3-2) emission 
at the 3-ct level, and use this to estimate a dense gas mass 
Of M d ense( H 2 )= qhcn LU < 4.5 x 10 10 M assuming 
qhcn = 10 (|Gao fc Solom on 2004). Within the uncertainties 
in QfHCN this is consistent with the gas mass derived from 
12 CO, indicating that HCN may be a promising route to 
constraining t he dense gas ma sses of high-redshift galaxies 
in the future |Gao et al. I2007I I. particularly if two or more 
HCN lines are measured. We also find that the HCN/FIR ra- 
tio is consistent with that of local ULIRGs and star-forming 
galaxies with no evidence for AGN contribution to the FIR 
luminosity. 

We use the atomic fine-structure and molecular line ra- 
tios to investigate the excitation conditions within the ISM, 
and find that they are similar to those found in local star- 
bursts and ULIRGs. Using a grid of PDR models, we show 
that the molecular gas in SMMJ2135 is exposed to a UV 
radiation field which is ~ 1000 x more intense than that 
of the Milky Way, and has a density of n ~ 10 4 cm -3 , 
both characteristic of the central regions of a typical star- 
burst galaxy. Thus, the inferred density and far-UV radia- 
tion field strength of the ISM appear similar to those seen in 
local ULIRGs consistent with the interpretation of SMGs as 
the high-redshift analogs of these merger-driven starbursts. 
However, the spatial extent and mass of the gas reservoir 
(and spatial extent of the star formation) is larger than typ- 
ically found in local ULIRGs. 

We find that the 12 CO emission shows multiple kine- 
matic components and we decompose the SLEDs into three 
kinematic components to investigate the excitation varia- 
tions within the system. We find tentative evidence of tem- 



perature variation between the kinematic components, al- 
though we defer a detailed discussion of the kinematics to a 
future paper (P. Cox et al. in prep). 

Using a number of approaches we have derived esti- 
mates for the effective size of the gas reservoir within this 
system. All of our estimates are larger than the size of the 
clumps detected in the 870-/im map from SMA, implying 
that the gas traces a more extended structure. Overall, our 
modelling of the 12 CO SLED supports the presence of two 
phases in the ISM, with one phase associated with ~solar 
metallicity star-forming clumps with > 100 pc sizes, em- 
bedded in a low metallicity (or high cosmic ray flux), ex- 
tended component with a characteristic size of ~ 1 kpc. 
However, even the "extended" component has a much higher 
gas density than typically observed in local starbursts such 
as M 82 (although a similar density to local ULIRGs such as 
Arp220), and is exposed to a UV radiatio n field which is ~ 
1000 x more intense than the Milky Way l|Gerin fc P hillips 
Il998ft . 

Assuming that the molecular gas is evenly distributed, 
we show that the integrated properties of the galaxy follow 
the Kennicutt-Schmidt relation. We also test where the indi- 
vidual star-forming regions would lie if they too followed the 
Kennicutt-Schmidt relation, and we show that on ~ 100 pc 
scales, the gas density would be ~ 30 x higher than that 
inferred from the LVG or PDR analysis. However, the prop- 
erties of the diffuse gas disk suggest a characteristic density 
n ~ 0.2-1 x 10 4 cm 2 , which is comparable to that inferred 
from the LVG and PDR models. These results are consis- 
tent with a model in which the sub-millimetre emission pre- 
dominantly traces the dense cores of the star-forming re- 
gions, but that these clumps are simply embedded in a much 
more extended and diffuse structure which dominates the 
gas emission. These results may also imply a break-down 
of the Kennicutt-Schmidt relation on scales of individual 
GMCs (e.g. < 100 pc), although clearly only high-resolution 
millimetre observations will accurately constrain their sizes. 

These results show that it is possible to unravel the 
complex ISM physics within starburst galaxies at high red- 
shift at a level of detail which, until recently has only been 
possible in the z < 0.1 Universe. Our observations pro- 
vide an unique window into the physics of star-formation 
at z ~ 2 with signal-to-noise that would otherwise require 
the increased light grasp of the next generation facilities. 
Indeed this source provides insight into the science that will 
be routinely possible once ALMA reaches full science oper- 
ations. 



ACKNOWLEDGMENTS 

ALRD acknowledges an STFC studentship, AMS gratefully 
acknowledges a Sir Norman Lockyer Royal Astronomical So- 
ciety Fellowship, IRS acknowledges support from STFC. We 
are very grateful to Padeli Papadopoulos for extensive disc- 
cusion and comments on this manuscript. We also thank Es- 
telle Bayet, Carsten Kramer, Rowin Meijerink, Serena Viti, 
Paul van der Werf and Chris Wilson for a number of use- 
ful discussions. We thank the anonymous referee for useful 
comments which added to the clarity and content of this 
paper. The observations in this paper were carried out with 
the IRAM 30 m and the Plateau de Bure Interferometer. 



© 0000 RAS, MNRAS 000, 000-000 



16 Danielson et al. 



IRAM is supported by INSU/CNRS (France), MPG (Ger- 
many) and IGN (Spain). We also thank the ESO Director for 
generously granting DDT observations with SHFI as part of 
program 283.A-5014 which is based on data acquired with 
the Atacama Pathfinder Experiment (APEX). APEX is a 
collaboration between the Max-Planck-Institut fur Radioas- 
tronomie, the European Southern Observatory, and the On- 
sala Space Observatory. The Submillimeter Array is a joint 
project between the Smithsonian Astrophysical Observatory 
and the Academia Sinica Institute of Astronomy and Astro- 
physics and is funded by the Smithsonian Institution and 
the Academia Sinica. Zpectrometer observations were taken 
as part of program 09A-040. Zpectrometer is supported by 
National Science Foundation grants AST-0503946 and AST- 
070865, and by the National Radio Astronomy Observatory. 
The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under coopera- 
tive agreement by Associated Universities, Inc. 



REFERENCES 

Aalto S., Booth R. S., Black J. H., Johansson L. E. B., 

1995, A&A, 300, 369 
Aravena M., Carilli C, Daddi E., Wagg J., Walter F., 

Riechers D., Dannerbauer H., Morrison G. E., Stern D., 

Krips M., 2010, ApJ, 718, 177 
Baker A. J., Tacconi L. J., Genzel R., Lehnert M. D., Lutz 

D., 2004, ApJ, 604, 125 
Bayet E., Gerin M., Phillips T. G., Contursi A., 2004, A&A, 

427, 45 

Bayet E., Gerin M., Phillips T. G., Contursi A., 2006, A&A, 
460, 467 

Bennett C. L., Fixsen D. J., Hinshaw G., Mather J. C, 
Moseley S. H., Wright E. L., Eplee, Jr. et al 1994, ApJ, 
434, 587 

Blain A. W., Smail I., Ivison R. J., Kneib J. -P., Frayer 

D. T., 2002, PhR, 369, 111 

Bolatto A. D., Jackson J. M., Kraemer K. E., Zhang X., 

2000, ApJL, 541, L17 
Bradford C. M., Nikola T., Stacey G. J., Bolatto A. D., 

Jackson J. M., Savage M. L., Davidson J. A., Higdon S. J., 

2003, ApJ, 586, 891 
Brown W. R., Kewley L. J., Geller M. J., 2008, AJ, 135, 

92 

Bryant P. M., Scovillc N. Z., 1996, ApJ, 457, 678 

Carilli C. L., Daddi E., Riechers D., Walter F., Weiss A., 

Dannerbauer H., Morrison G. E., Wagg, J et al. 2010, 

ApJ, 714, 1407 
Coppin K. E. K., Swinbank A. M., Neri R., Cox P., Smail 

I., Ellis R. S., Geach J. E., Siana et al. 2007, ApJ, 665, 

936 

Dave R., Finlator K., Oppenheimer B. D., Fardal M., Katz 
N., Keres D., Weinberg D. H., 2010, MNRAS, 404, 1355 

Davies R. I., Sternberg A., Lehnert M., Tacconi- Carman 
L. E., 2003, ApJ, 597, 907 

Devereux N., Taniguchi Y., Sanders D. B., Nakai N., Young 
J. S., 1994, AJ, 107, 2006 

Downes D., Solomon P. M., 1998, ApJ, 507, 615 

Flower D. R., 2001, MNRAS, 328, 147 

Flower D. R., Le Bourlot J., Pineau Des Forets G., Roueff 

E. , 1994, A&A, 282, 225 



Frayer D. T., Ivison R. J., Scoville N. Z., Evans A. S., Yun 

M. S., Smail I., Barger A. J., Blain A. W., Kneib J.-R, 

1999, ApJL, 514, L13 
Frayer D. T., Ivison R. J., Scoville N. Z., Yun M., Evans 

A. S., Smail I., Blain A. W., Kneib J., 1998, ApJL, 506, 

L7 

Frerking M. A., Keene J., Blake G. A., Phillips T. G., 1989, 
ApJ, 344, 311 

Gao L., Springel V., White S. D. M., 2005, MNRAS, 363, 
L66 

Gao Y., Carilli C. L., Solomon P. M., Vanden Bout P. A., 

2007, ApJL, 660, L93 
Gao Y., Solomon P. M., 2004, ApJ, 606, 271 
Genzel R., Tacconi L. J., Gracia-Carpio J., Sternberg A., 

Cooper M. C, Shapiro K., Bolatto A., Bouche N., Bour- 

naud F., Burkert, A., et al. 2010, ArXiv e-prints 
Gerin M., Phillips T. G., 1998, ApJL, 509, L17 
Gerin M., Phillips T. G., 2000, ApJ, 537, 644 
Gonzalez- Alfonso E., Fischer J., Isaak K., Rykala A., Savini 

G., Spaans M., van der Werf P., Meijerink, R. et al 2010, 

ArXiv e-prints 

Greve T. R., Bertoldi F., Smail I., Neri R., Chapman S. C, 
Blain A. W., Ivison R. J., Genzel et al. 2005, MNRAS, 
359, 1165 

Greve T. R., Ivison R. J., Papadopoulos P. P., 2003, ApJ, 
599, 839 

Greve T. R., Papadopoulos P. P., Gao Y., Radford S. J. E., 
2009, ApJ, 692, 1432 

Guesten R., Serabyn E., Kasemann C, Schinckel A., 
Schneider G., Schulz A., Young K., 1993, ApJ, 402, 537 

Hailey-Dunsheath S., Nikola T., Stacey G. J., Oberst T. E., 
Parshley S. C, Benford D. J., Staguhn J. G., Tucker C. E., 
2010a, ApJL, 714, L162 

Hailey-Dunsheath S., Nikola T., Stacey G. J., Oberst T. E., 
Parshley S. C, Benford D. J., Staguhn J. G., Tucker C. E., 
2010b, ApJL, 714, L162 

Harris A. I., Baker A. J., Jewell P. R., Rauch K. P., 
Zonak S. G., O'Neil K., Shelton A. L., Norrod, R. 
et al., 2007, in A. J. Baker, J. Glenn, A. I. Harris, 
J. G. Mangum, & M. S. Yun ed., From Z-Machines 
to ALMA: (Sub)Millimeter Spectroscopy of Galaxies 
Vol. 375 of ASPC, The Zpectrometer: an Ultra- Wideband 
Spectrometer for the Green Bank Telescope, p. 82 

Harris A. I., Baker A. J., Zonak S. G., Sharon C. E., Genzel 
R., Rauch K., Watts G., Creager R., 2010, ArXiv e-prints 

Harrison A., Puxley P., Russell A., Brand P., 1995, MN- 
RAS, 277, 413 

Hocuk S., Spaans M., 2010, A&A, 510, A110+ 

Israel F. P., Baas F., 2002, A&A, 383, 82 

Israel F. P., Maloney P. R., Geis N., Herrmann F., Madden 
S. C, Poglitsch A., Stacey G. J., 1996, ApJ, 465, 738 

Israel F. P., White G. J., Baas F., 1995, A&A, 302, 343 

Ivison R. J., Smail I., Papadopoulos P. P., Wold I., Richard 
J., Swinbank A. M., Kneib J., Owen F. N., 2010, MNRAS, 
404, 198 

Ivison R. J., Swinbank A. M., Swinyard B., Smail I., Pear- 
son C. P., Rigopoulou D., Polehampton E., Baluteau, J.-P. 
et al. 2010, A&A, 518, L35 

Juneau S., Narayanan D. T., Moustakas J., Shirley Y. L., 
Bussmann R. S., Kennicutt R. C, Vanden Bout P. A., 
2009, ApJ, 707, 1217 

Kaufman M. J., Wolfire M. G., Hollenbach D. J., Luhman 



© 0000 RAS, MNRAS 000, 000-000 



The ISM in a z ~ 2.3 Star-Forming Galaxy 17 



M. L., 1999, ApJ, 527, 795 
Kennicutt Jr. R. C, 1998, ApJ, 498, 541 
Kramer C, Cubick M., Rollig M., Sun K., Yonekura Y., 

Aravena M., Bensch F., Bertoldi, F. et al. 2008, A&A, 

477, 547 

Kramer C, Jakob H., Mookerjea B., Schneider N., Briill 
M., Stutzki J., 2004, A&A, 424, 887 

Krumholz M. R., Tan J. C, 2007, ApJ, 654, 304 

Maiolino R., Caselli P., Nagao T., Walmsley M., De Breuck 
C, Meneghetti M., 2009, A&A, 500, LI 

Malhotra S., Kaufman M. J., Hollenbach D., Helou G., Ru- 
bin R. H., Brauher J., Dale D., Lu, N. Y. et al 2001, ApJ, 
561, 766 

Mao R. Q., Henkel C, Schulz A., Zielinsky M., Mauers- 
berger R., Storzer H., Wilson T. L., Gensheimer P., 2000, 
A&A, 358, 433 
Meijerink R., Spaans M., Israel F. P., 2007, A&A, 461, 793 
Negishi T., Onaka T., Chan K., Roellig T. L., 2001, A&A, 
375, 566 

Onodera S., Kuno N., Tosaki T., Kohno K., Nakanishi K., 
Sawada T., Muraoka K., Komugi, S. et al. 2010, ArXiv 
e-prints 

Panuzzo P., Rangwala N., Rykala A., Isaak K. G., Glenn 
J., Wilson C. D., Auld R., Baes, M. et al. 2010, ArXiv 
e-prints 

Papadopoulos P. P., 2007, ApJ, 656, 792 
Papadopoulos P. P., Greve T. R., 2004, ApJL, 615, L29 
Papadopoulos P. P., Isaak K. G., van der Werf P. P., 2007, 
ApJ, 668, 815 

Papadopoulos P. P., Thi W., Viti S., 2004, MNRAS, 351, 
147 

Papadopoulos P. P., van der Werf P., Isaak K., Xilouris 
E. M., 2010, ApJ, 715, 775 

Pety J., 2005, in F. Casoli, T. Contini, J. M. Hameury, 
& L. Pagani ed., SF2A-2005: Semaine de l'Astrophysique 
Francaise Successes of and Challenges to GILD AS, a 
State-of-the-Art Radioastronomy Toolkit, p. 721 

Pineda J. L., Bensch F., 2007, A&A, 470, 615 

Sakamoto K., Scoville N. Z., Yun M. S., Crosas M., Genzel 
R., Tacconi L. J., 1999, ApJ, 514, 68 

Scoville N., 2004, in S. Aalto, S. Huttemeister, & A. Pedlar 
ed., The Neutral ISM in Starburst Galaxies Vol. 320 of 
ASPC, Starburst and AGN Connections, pp 253-+ 

Scoville N. Z., Yun M. S., Bryant P. M., 1997, ApJ, 484, 
702 

Smail I., Ivison R. J., Blain A. W., Kneib J.-R, 2002, MN- 
RAS, 331, 495 

Socrates A., Davis S. W., Ramirez-Ruiz E., 2008, ApJ, 687, 
202 

Solomon P. M., Vanden Bout P. A., 2005, ARAA, 43, 677 

Spaans M., 1996, A&A, 307, 271 

Spaans M., Meijerink R., 2008, ApJL, 678, L5 

Spergel D. N., Verde L., Peiris H. V., Komatsu E., Nolta 

M. R., Bennett C. L., Halpern M., Hinshaw, G., et al. 

2003, APJS, 148, 175 
Stacey G. J., Geis N., Genzel R., Lugten J. B., Poglitsch 

A., Sternberg A., Townes C. H., 1991, ApJ, 373, 423 
Stark A. A., Bolatto A. D., Chamberlin R. A., Lane A. P., 

Bania T. M., Jackson J. M., Lo K., 1997, ApJL, 480, L59 
Stark D. P., Swinbank A. M., Ellis R. S., Dye S., Smail 

I. R., Richard J., 2008, Nature, 455, 775 
Stutzki J., Graf U. U., Haas S., Honingh C. E., Hottgenroth 



D., Jacobs K., Schieder R., Simon, R. et al 1997, ApJL, 
477, L33 

Swinbank A. M., Smail I., Longmore S., Harris A. I., Baker 

A. J., De Breuck C, Richard J., Edge, A. C. et al 2010, 

Nature, 464, 733 
Tacconi L. J., Genzel R., Smail I., Neri R., Chapman S. C, 

Ivison R. J., Blain A., Cox et al. 2008, ApJ, 680, 246 
Tacconi L. J., Neri R., Chapman S. C, Genzel R., Smail 

I., Ivison R. J., Bertoldi F., Blain et al. 2006, ApJ, 640, 

228 

Teyssier R., Chapon D., Bournaud F., 2010, ArXiv e-prints 
Thompson T. A., Quataert E., Murray N., 2005, ApJ, 630, 
167 

van der Werf P. P., Isaak K. G., Meijerink R., Spaans M., 
Rykala A., Fulton T., Loenen A. F., Walter, F 2010, ArXiv 
e-prints 

Vassilev V., Meledin D., Lapkin I., Belitsky V., Nystrom 
O., Henke D., Pavolotsky A., Monje, R. et al. 2008, A&A, 
490, 1157 

Walter F., Riechers D., Cox P., Neri R., Carilli C, Bertoldi 
F., Weiss A., Maiolino R., 2009, Nature, 457, 699 

Ward J. S., Zmuidzinas J., Harris A. I., Isaak K. G., 2003, 
ApJ, 587, 171 

Weifi A., Downes D., Henkel C, Walter F., 2005, A&A, 
429, L25 

Weifi A., Downes D., Neri R., Walter F., Henkel C, Wilner 
D. J., Wagg J., Wiklind T., 2007, A&A, 467, 955 

Weifi A., Downes D., Walter F., Henkel C, 2005, A&A, 
440, L45 

Weifi A., Walter F., Scoville N. Z., 2005, A&A, 438, 533 
White G. J., Ellison B., Claude S., Dent W. R. F., Math- 

eson D. N., 1994, A&A, 284, L23 
Wild W., Harris A. I., Eckart A., Genzel R., Graf U. U., 

Jackson J. M., Russell A. P. G., Stutzki J., 1992, A&A, 

265, 447 

Young J. S., Scoville N. Z., 1984, ApJ, 287, 153 
Yuan T., Kewley L. J., 2009, ApJL, 699, L161 
Zmuidzinas J., Betz A. L., Boreiko R. T., Goldhaber D. M., 
1988, ApJ, 335, 774 



© 0000 RAS, MNRAS 000, 000-000 



